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FOREWORD 

In the summer of 1983, Mr. C.L. Ball arranged for the authors to 
spend three weeks in residence at NASA-Lewis Research Center, chiefly to 
consider together how to combine recent theories of surge and rotating 
stall. The present report describes the results of that effort. The 
authors would like to acknowledge Mr. Ball's assistance, and also the 
contribution of Mr. R. Chue of MIT who provided certain supporting 
calculations described in Chapter 3 of this report. 
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1 . INTRODUCTION 
Background 

Problems of compressor instability have been of concern to aircraft 
engine designers for a number of years, and the provision of sufficient 
stall margin is an important consideration in any compressor design. In 
addition to problems associated with the inception of stall, however, 
another aspect of this general topic has been of increasing import in 
recent years. This is the behavior of the compression system subsequent to 
the onset of stall, i.e., subsequent to the onset of compressor or 
compression system instability. 

The reason for this increased interest in post-stall behavior has to 
do with the phenomenon of non-recoverable stall (or stall "stagnation"). 

In this condition (described in [1]) the hysteresis associated with 
compressor operation in rotating stall is large enough that, once the 
engine has encountered rotating stall, it cannot recover to a condition of 
unstalled operation. The only remedy then is to greatly decrease the 
engine speed to decrease the hysteresis, or, in some cases, to shut the 
engine down and restart it. 

A key task associated with this problem is to predict post-stall 
behavior. This knowledge is essential for rational design of 
stagnation-resistant compressors. This is true not only for quasi-steady 
operation but also for the unsteady features of the compression system 
response. In the first instance, one may be able to consider the 
compressor in isolation, without coupling it to the system. This is 
essentially the problem treated by Day, Greitzer, and Cumpsty [2] 
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(semi-empirically ) , and, more recently, by Moore [3] in a rather more 
fundamental manner. To grapple with the actual situation in more 
generality, however, one must deal with the unsteady behavior of the 
compressor during the mass flow oscillations that characterize post-stall 
transients. This involves a description of the possible formation (growth) 
in the compressor of a rotating stall cell during this period, while at the 
same time, the compressor interacts with other components of the system. 

It is to be stressed from the outset that this task cannot today be 
carried out by exact or numerical solution of the basic partial 
differential equations of fluid motion. The prediction of even steady- 
state rotating stall performance has not yet been carried out on this 
level. However, many of the salient features of rotating stall have been 
predicted using certain simplifying assumptions, in a theory by Moore [3]. 
We do not hope to unravel all of the issues connected with this flow 
regime, but it is nevertheless felt that a useful model for transient 
compressor operation can be formulated in a similar way. This will not 
only be helpful for general physical insight, but should give guidance to 
experimental studies, and also illustrate some of the overall effects of 
different design parameters. This report, which essentially covers the 
work done by the authors during the period 7/25/83 to 8/12/83 while both 
were in residence at NASA Lewis Research Center, describes an approximate 
theory of this type. 

Motivation: Time Scales in Transient Compressor Operation 

Before plunging into the details of the analysis, it seems appropriate 
to provide a motivating discussion of the physical situation. This is 
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perhaps especially true in view of the large amount of work that has been 
done on dynamic behavior of compressors and compression systems. In this 
section, therefore, we will present some basic physical arguments that 
underlie an examination of the problem under consideration; that is, 
combined surge and rotating stall. In particular, we will show that there 
appears to be a strong reason for an examination of the transient behavior 
of the compressor as it passes in and out of rotating stall, since its 
behavior can affect the overall system response. 

One way to do this is to consider the time scales that are associated 
with the phenomena that occur during the transients of interest. The first 
of these is the rotating stall formation time. This is the time needed 
for the flow in the compressor to change from a nominally uniform 
(unstalled) to a severely asymmetric (rotating stall) state. This change 
can be associated with the shedding of the bound circulation of a 
significant portion of the blading in the compressor (say, one-third to 
one-half) and the subsequent convection of this shed vorticity downstream 
on the order of the circumference. The time taken for this to occur might 
thus be on the order of the flow time for this distance, i.e., ttD/C x 
where D is the diameter of the compressor and C x is some representative 
axial velocity. 

The above argument suggests the use of itD/C x as the proper rotating 
stall formation time. In the past, however, this time has been quoted in 
terms of rotor revolutions, i.e., in terms of ttD/U where U is the blade 
speed. While perhaps not as physically grounded as ttD/C x , available data 
is in terms of this latter parameter. Also, once a given compressor is 
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specified, the axial velocity and the rotor speed are of similar magnitudes, 
linked by geometry. Thus we will defer to this convention and specify 

Rotating stall formation time scale = t = ttD/U (1.1) 

K. 

as the time scale for rotating stall cell formation. Available data 
indicates that the characteristic time is roughly two rotor revolutions 
and that the time to "fully develop" is even longer, being several 
characteristic times. 

Another relevant time scale is the "flow change time." This is the 
characteristic time for changes in the overall (annulus averaged) mass 
flow through the compressor. It is set by a balance between the pressure 
forces that are available to accelerate the flow in the compressor duct 
during a stall transient and the inertial forces that arise from this 
acceleration. 

If the overall compressor pressure rise is AP and the pre-stall axial 
velocity is denoted by C x then this balance of pressure and inertial 
force yields 
L C 

P = KAP (1.2) 

c 

where KAP is a number of order unity. In this expression the quantity 
L c is the effective length of the compressor duct and x c is the flow 
change time. We can write this in terms of non-dimensional compressor 
parameters $( = C x /U) and \l>( = AP/pU^) to yield an expression for the 
flow change time scale as 
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L L 

C 6 „ C 

Flow change time scale = t = — rr~ - Z — 
^ c U Kt|j U 


( 1 . 3 ) 


In this expression Z is a non-dimensional number that is a function of the 
compressor parameters* Z is expected to be of order magnitude unity* 

A third time scale is the "Helmholtz resonator time." This is the 
characteristic time for small oscillations in a compression system, 
considered as a Helmholtz resonator* Using the formula for the Helmholtz 


frequency, this is 


Helmholtz resonator time scale ~ t = — = 


i V L 
1 p c 

a A 

s c 


where Vp is the system plenum volume and A c is the through-flow area of 
the compressor. 

For low speed machines or for compressors of few stages, the large 
amplitude (surge) system oscillation time scale may also be of this order 
so that it and the flow change time will be roughly the same* However, 
for multistage machines, especially at high speed, the transients that 
occur subsequent to stall are not at all sinusoidal, but are rather of the 
nature of relaxation oscillations. Hence the flow change time scale and 
the Helmholtz resonator time scale can be quite different. (Note that for 
the relaxation type of transient, the plenum pressure change during the 
transient is small because the transient is so rapid. Thus the actual 
plenum volume should not enter into the definition of this time scale.) 

There is also a fourth time scale which we also include for 
completeness. This is the "plenum emptying time." It is only 
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appropriate to the relaxation type of oscillation, as it is set by a 
balance between resistive and capacitive elements, i.e., between flow 
resistance in the throttle and the pressure difference (plenum to 
atmosphere) due to the storage of potential energy of compression of the 
gas in the plenum. It can be much longer than the flow change time and is 
the principle part of the period of surge cycles experienced by multistage 
compressors near design speed. Explicitly this time is given by: 

V V 

Plenum emptying time = t = — = E(— ) — M c \ 

e <j>aaA a a A 1 * ' 

s s c s sc 


It is useful to put the foregoing times into non-dimensional ratios. 
These will contain not only the non-dimensional numbers Z and E introduced 
above, but two other non-dimensional parameters as well. These are: 


= . - u = 

2u L 2a »A L 
He s c c 


which has been described in detail in [3] and [4], and 
2L 

c 


( 1 . 6 ) 


(1.7) 


the ratio of compressor effective length to radius. 

The non-dimensional ratios of the various time scales are thus: 


Helmholtz resonator time T H 

Rotating stall formation time t 


Flow change time 

Rotating stall formation time 2ir 


Plenum emptying time 

Rotating stall formation time 



IT 


C 


( 1 . 8 ) 


(1.9) 


( 1 . 10 ) 
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If one substitutes representative numbers into the above non- 
dimensional ratios, it becomes evident that the first two at least are 
roughly of magnitude unity. In other words, in practical circumstances 
the rotating stall cell formation time is of the same order as the mass 
flow change time. This is the central point of this discussion concerning 
time scales. 

The implication of this rough parity (between the time scale of 
rotating stall cell formation and the time for a significant change in the 
overall mass flow during a transient) is that rotating stall will not have 
time to fully form during this type of non-steady flow. This has 
significant implications in turn for the modelling of post-stall 
transients, since, during the transient, the compressor is not performing 
in a quasi-steady manner. Thus, information obtained from a steady-state 
compressor test with the machine operating in rotating stall may not be at 
all representative of the performance during a rapid post-stall transient. 
If so, it is imperative to address the question of compressor behavior 
during these transients. Specifically, the theory developed herein must 
include an assessment of how the growth or decay of rotating stall affects 
the instantaneous compressor pumping characteristic, and hence overall 
system behavior. 

Before leaving the discussion of time scales, there is one further 
point that should be mentioned. The principle aspect of compressor 
unsteady response that we are examining is that associated with the growth 
of an asymmetric flow pattern which has a length scale on the order of the 
compressor diameter. The time scale is thus very much longer than that 
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associated with the unsteady response of a single blade passage. This 
latter might be on the order of the through-flow velocity, W, divided by 
b, the blade chord, and hence the ratio of the blade response time to the 
rotating stall formation time is of order 

Blade through flow time ~ b/W _ b U 

Rotating stall formation time - irD/U 7 d W <<: 1 (1.11 

Thus, the precise modelling of the unsteady performance of the blades, per 
se, would not be expected to be a key ingredient of a theory of stall-cell 
growth. 

Overall Format of the Report 

In the next section we discuss the physical assumptions inherent in 
the compressor and compression system models that we will use. We then 
briefly describe separately the two "components" of the overall problem, 
rotating stall and surge. An analysis is first presented of rotating 
stall (where the annulus averaged flow is constant in time, but the flow 
is non-uniform around the circumference). The pure surge mode (in which 
the annulus averaged mass flow varies with time, but the flow is uniform 
around the circumference) is then examined. Approximate solutions are 
developed for both of these modes. 

Following these two sections, which are to some extent reviews of 
previous work that the present analysis builds on, the combined problem is 
introduced. Using an approximate (Galerkin) procedure, a set of coupled 
nonlinear equations is derived which describes the time varying state of 
the system. The nature of the equations is discussed, and then certain 
numerical solutions are described which show the effects of certain para- 
meters of interest. Finally, some comments are made about future avenues 
of research, both theoretical and experimental, that are thought to be 
needed for this general problem. 
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2. FLUID DYNAMIC MODEL 

The basic compression system to be analyzed is shown in Figure 2.1. 
The coordinate system that will be used, and the directions of the 
different velocity components, as well as all the different axial stations 
that we will employ, are also indicated. 

The system consists of a compressor operating in a duct and 
discharging to a downstream plenum. The plenum dimensions are large 
compared to those of the compressor duct so that the velocities and fluid 
accelerations in the plenum can be considered negligible. Hence the 
pressure in the plenum can be considered to be uniform spatially, although 
varying in time. The flow through the system is controlled by a throttle 
at the plenum exit. 

The systems that we consider will have overall pressure rises, 
atmosphere to plenum, that are small compared to the ambient level. If, 
also, we assume Mach numbers to be small, and oscillations to have 
frequencies well below those of acoustic resonance, then we can adopt an 
incompressible flow description of compressor behavior (although, of 
course, compressibility is important in setting system dynamics). The 
basic compression system model, as discussed so far and shown in Figure 
2.1, has been used by many investigators. Where we depart from previous 
practice in analysis of compression system transients, however, is in the 
coupling of an analysis of the two-dimensional unsteady flow in the 
compressor with the lumped-parameter system model. We therefore should 
also briefly describe our representation of the compressor. 

The basic compressor model is similar to that presented in [3]. The 
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compressor is considered to have a high hub/tip radius ratio, so that a 

two dimensional description can be used. The compressor duct is modeled as 

having an inlet section (contraction) followed by a constant-area section 

upstream of the compressor. The annulus immediately downstream of the 

compressor also has constant area. These restrictions simplify the 

analysis, but could be removed in a more general treatment. 

The following nomenclature will be used (see Figure 2.1): all 

distances will be non-dimensionalized by the mean radius. The 

non-dimensional circumferential and axial coordinates and time are: 

Circumferential coordinate: e - circumferential distance _ _ 

j < 2 -" 

Axial coordinate: n = axial distance/(R) (2.2) 

Time: £ = Ut/R (2.3) 

As discussed in [3], the basic model we adopt for the unsteady 
performance of a compressor blade row is that the pressure rise across a 
single row is given by 


AP 



2 


F(<t>) - 


dcfr 

dt 


(2.4) 


where <{>(= C x /U) is the axial velocity coefficient at the compressor. The 
quantity F(<j)) is the axisymmetric performance characteristic, and t can be 
viewed as a "time constant" associated with the internal lags in the 
compressor. In [3], it was shown that a reasonable value for t is roughly 
t ~ 2b x /C x . 

The unsteadiness in the flow through a stator passage reflects the 
accelerations associated with transient effects. For a rotor, however, 
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there is another source of the unsteadiness seen by a blade row due to the 
rotor blades moving (with velocity U) through a circumferentially 
non-uniform flow* Therefore, 


d<|> \ _ 3 _ U 3 

■ stator ^ 


m 

' dt / rotor 


= 

3 1 


U 3j> __ U / 3<j) 
+ R 30 _ R \3? 



(2.5) 

( 2 . 6 ) 


We note that in the present situation, because of the possibility of growth 
or decay of velocity non-uniformities, there is not necessarily Galilean 
equivalence between the changes in time and the spatial derivatives, as 
there would be for a pure traveling wave. 

Proceeding to the multistage situation, we have to give the pressure 
rise across an N-stage compressor (not including the inlet and exit guide 
vanes which we will discuss subsequently): 

^ (* S ♦ S) 

pu 

Rather than introduce more of the detailed features of the model at 
this point, we will next examine the two special cases that were mentioned 
previously, with discussion of further aspects of the model being intro- 
duced as needed. Since the development of the pure surge and rotating 
stall situations have appeared elsewhere, however, we will not present the 
derivation ab initio, but rather refer the reader to [3], [4], [5], and [6] 
for fuller accounts. 
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3. ROTATING STALL ANALYSIS 

Governing Equation for Rotating Stall (Constant annulus averaged mass flow) 
We begin with the equation governing a disturbance (g) of axial flow 
coefficient ($) 

<t> ® $ + g (3.1 ) 

where $ is the steady mean flow coefficient (or nondimensional axial 
velocity) and g is the nonaxi symmetric disturbance. This nonaxisymmetric 
part of the axial velocity will be a function of 0*, where e* = 0 - f £ . 

That is, 0* is measured in a coordinate system traveling with the stall 
cell, which propagates around the circumference with a nondimensional, 
constant speed f. (The physical speed in the laboratory coordinate 
system is fU. ) 

Equation (1) of Part III of [3] is 

- mfh + ¥ - \|> (<|>) = 0 (3.2) 

de c 

where h is the nondimensional circumferential velocity distrubance. Like 
g, h is also a function of 6*. The quantity m is a parameter which 
reflects inlet and outlet flow inertia effects. If the exit flow is a 
sudden expansion, and entrance flow is potential, then m = 1. If the exit 
is a straight duct, we know only that m = 2 for small disturbances. We 
assume that m is constant for large disturbances as well, with m = 2 being 
appropriate for long diffusers. The parameter X depends on f and a, the 
latter being a dimensionless expression of the time lag (t) of pressure 
rise for a single blade passage, as described in [3]. 
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X 


i <i - f > 


(3.3) 


a = R/(UNt) (3.4) 

The quantity a is assumed known while f (and hence X) is to be found as an 


eigenvalue of the solution. 

It is emphasized in [3] that although the compressor characteristic 
function i|>^.(<()) is assumed known, it must be identified as an "axi symmetric 
characteristic," applicable when there is no rotating stall. The actual 
steady pressure-rise coefficient (Y) produced when the compressor is in 
rotating stall* differs from ip c ( <J> ) because there is a nonlinear effect of 
rotating stall itself. That is, even when <j> = $ (the average flow 
coefficient), Y and \|> c differ. We denote that difference by 6: 

5 = Y - i|> {$) (3.5) 

c 


so the last two terms of Eq. (3.2) are 

-[ib (d>) - i|j ($)] + 6 = -G(g) + 6 (3.6) 

c c 

and the quantity in brackets, depending only on <(>-$, or g, we write as 
G(g) . These relationships are illustrated on Fig. 3.1. Equation (3.2) 
may now be written 

X% - mfh - G(g) +6=0 (3.7) 

de 

It is necessary to specify how h(0*) and g(0*) are related before 
proceeding to solve Eq. (3.7). 


*This would be the coefficient measured during a steady-state compressor 
test. 
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The "dh/de* = -g" Approximation 

The relation between h and g is discussed in Appendix A of [4], in 


general, complete flow field solutions in the entrance and exit regions are 
needed to relate h and g, and no simple exact relationship applies except 
in special cases. In [3] and [4] arguments are given, however, to support 
the adoption of the approximate relation 



(3.8) 


Some further discussion of this point is also presented in Appendix A of 
the present report. The general conclusion from these discussions is that 
the approximation can be adopted for the nonlinear stall cell calculation 
with some confidence. We will therefore do this in the remainder of the 
report. 

If this is done, the set of equations that describe pure rotating 
stall can be written by introducing Eq. (3.8) into the derivative of Eq. 
(3.7) to yield 


d g dG dg 

X - — — - + mfg = 0 

de*2 d 9 de* 


(3.9) 


The axial velocity perturbation, g, is also required to be periodic with 
period 2ir. 

g(9* + 2ir) = g( 6* ) 

The cyclic integral of h is required to vanish, in order that the 
transverse velocity disturbance introduces no circulation. 
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2tt 

hde = 0 

0 


(3.11) 


The condition given in Eq. (3.11) is strictly valid only when the entrance 
flow is potential. Further research will be needed to explore conditions 
alternative to Eq. (3.11). In view of Eqs. (3.10) and (3.11), integration 
of Eq. (3.9) over a full cycle yields the performance effect of rotating 
stall: 


6 


_ 1 _ 

2t r 


2ir * 
G (g )d0* 

0 


(3.12) 


Equations (3.9), (3.10), and (3.12) will represent the basic rotating- 
stall problem: a given throttle setting defines the origin of g and allows 

the given axisymmetric characteristic to be expressed (relative to the 
throttle point) as the function G(g). Then, the eigensolution g(6*) is 
found by integrating Eq.(3.9) with f (and hence X) found as an eigenvalue 
to fix the period at 2 tt . The nonlinear performance effect then follows 
from Eq. (3.12) . 

Numerical Solution for a Cubic Characteristic 

In Refs. [7] and [8] it is argued that the axisymmetric characteristic 
is typically a smooth S-shaped curve. (In those references, the axisym- 
metric curve is estimated for a three-stage compressor using a procedure 
based on corrections to transient data.) If this idea is correct, a 
physically realistic function would be a simple cubic, as that shown in 
Fig. 3.2. We might therefore expect that the in-rotating-stall character- 
istic calculated using the cubic will be a good representation of that 
calculated using an exact axisymmetric diagram. 
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The formula for the particular curve shown in Fig. 3.2 is 

* c (<|)) = + H [1 + f ( i " 1} “ ^ " 1 ) 3 J (3.13) 

o 


We locate the throttle-line intersection at a point B to the right of the 
midpoint of the diagram. Since we chose to measure <|> from the minimum of 
i|» c (representing the shut-off condition in this case), and since G and g 
are both measured from the throttle-line intersection, 


G = i|> U) - t (W + B) 
c c 


(3.14) 


g = <J> - (w + b) 

Substituting into Eq. (3.21) yields 


G ( g) - [3 (1 - \)g - 3^-g 2 

W w 



] 


(3.15) 


(3.16) 


which we will use hereafter to represent the damping function G(g). 

In Refs. [7] or [8], a numerical procedure is presented for solving 
Eqs. (3.9) to (3.12), and as part of the present study, R. Chue of MIT 
carried out calculations for a particular cubic G(g). The value used for m 
was 1.75. H and W were chosen so that the cubic would have the same peak 
and valley points as the actual curve; their values were 0.18 and 0.25 
respectively. Rather than Eq. (3.3), a somewhat more elaborate relation 
between X and f was used, in which guide-vane effects were included, 
following Ref. [8]. Numerically, the relation was, in effect, 

X = 4.34(0.378 - f) 

Figure 3.3 shows how f increases when the throttle is opened (i.e., as 
B is increased), as predicted in [3]. Figs. 3.4a-f show how the axial 
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velocity disturbance g(6*) evolves as the throttle is opened. When B = 0, 
the symmetry features shown in Fig. 3.4a must occur, and we note that g(e*) 
swings between extremes which lie on the unstalled and reverse-flow legs of 
the characteristic. As B increases, these extrema remain nearly the same, 
but more time is spent unstalled, so to speak. There is no solution beyond 
B = 0.25, which corresponds to the peak point of the curve. 

The most interesting result is shown in Fig. 3.5. There, the cubic 
axisymmetric characteristic is shown with the calculated 6 superposed to 
form the predicted in-stall, pressure-rise variation with throttle setting 
($). The graph dramatically illustrates the idea of Ref. [7], that the 
usually-observed sudden break in performance at the stall point, and 
subsequent drop into M deep stall" as the throttle is closed, is the result 
of being in rotating stall, rather than a direct effect of the inherent 
axisymmetric characteristic of the compressor. The latter might, in fact, 
be quite smooth and gradual, with no break at all. 

Approximate Solution Using Galerkin Procedure 

Because we shall be concerned with the analysis of complicated oscill- 
ations of which rotating stall will be only a special case, it is useful, 
at this initial stage of the investigation, to devise a simple method of 
analytically representing rotating stall alone, with the idea of then using 
it later for more general post-stall disturbances. If the damping function 
G (g ) is smooth, the Galerkin method of nonlinear mechanics (Ref. [9]) is 
perhaps appropriate. (In this connection we note that this technique has 
been used quite successfully for the van der Pol problem, which is just the 


symmetric version of the present situation with B = 0.) 
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In general, in this procedure, the solution to the differential 
equation is represented by a suitable sequence of basic functions. If the 
functions are well chosen, and enough terms in the sequence are taken, the 
"true" solution can be represented very accurately. Fourier series or 
spectral methods are special cases of this procedure. 

We are seeking, for present purposes, the simplest possible wave 
representation capable of describing the stall phenomena. We have thus 
taken the somewhat drastic step of truncating the sequence after one 
harmonic term; in effect, a single term Galerkin procedure. The form 
chosen for g is therefore 

g cc sin 0* 

For the present exploratory purposes, this approximation will prove to 
be reasonable. However, it is very important to be clear about the limita- 
tions of this approximation. The single-term Galerkin procedure is most 
accurate for weakly nonlinear systems. The present set of equations is 
strongly nonlinear. Thus, although it will turn out that the approximation 
for g, the axial velocity disturbance, is reasonable, as are the pressure 
rise curves, which are averaged functions of g, the method does not give a 
good description of the derivatives of g. This is to be expected since the 
oscillations represented by Eq. (3.9) with small X would be of the 
relaxation type, and it would need many more than one harmonic to describe 
the shape accurately. 

In this study, we wish to emphasize the amplitudes of disturbances and 
integrated effects on performance. Thus, the information that is of most 
interest at present is just that information which is least sensitive to 
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the higher harmonics of the stall-cell description. Hence, for our 
purposes, the single-term approximation should be adequate, and should lead 
to a compact description of the overall stall behavior. Appendix B gives 
some results obtained using a two-term sequence of functions, as well as 
some additional discussion on the effects of incorporating more of the 
complete set of the basic functions. 

We now take g to be of the specific form 

g = A*W sin 6* . (3.17) 

Substituting into Eqs. (3.7, 3.16) yields 

i He*) = (-X + mf)A*sine* [(1-B 2 /W 2 )A* - ( 1/4)A* 3 ]cos0* (3.18) 

If the equation were satisfied exactly, the residual L(e*) would be 0. It 
is not, and L(e*) represents the error involved. In keeping with the 
Galerkin procedure, we require that the sine* and cose* moments of L(e*) 
vanish. This yields expressions for stall-cell speed (f) and wave ampli- 
tude (A*): 

X = mf (3.19) 

and 

A* = 2(1 - 0 2 /W 2 ) 1/2 (3.20) 

The form of g is thus 

g = [2W(1 - B 2 /W 2 ) 1/2 ] Sine* (3.21) 

Note that in the one-harmonic solution, there is no question as to the 
dh/de 1 = -g approximation, because only a single term is used. 

Equation (3.21) can be substituted into Eqs. (3.12, 3.17) to give the 
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result for <5 


6 = -3H( B/W) (1 - 8 2 /W 2 ) 


(3.22) 


Hence, according to Eq. (3.5) the rotating stall performance is 

¥ = + H(1 - 4 — + — — ) 

v c 2 W 2 3 

o W 


(3.23) 


where we recall that the parameter 8 represents the effect of the throttle 
setting, as shown on Fig. 3.2. This last formula can be compared with the 
axisymmetric characteristic (Eq. (3.13)): 


* c = * c + H(1 + 


3 B 1£1 
2 W " 2 .r 
w 


(3.24) 


Comparison with numerical solution . Equation (3.24) has already 
been indicated on Figs. 3.2 and 3.5. On Fig. 3.5 we also show Eq. (3.23), 
which may be compared with the exact result. Although the two curves do 
not agree numerically, it is remarkable that the qualitative features of 
the performance effect are well predicted by the simple, first-harmonic 
theory, especially the steep rise near 8 = W, and the equally steep drop 
near 8 = -W. 

Figure 3.6 shows the amplitude of the disturbance as a function of 
8/W. It vanishes at the two extremes of the diagram, 8 = ±W, and has its 
greatest value (2W) at the midpoint, 8=0. Although it is satisfying that 
"recovery" (in the sense of a limiting throttle setting) at 6 = W is very 
nearly correct, it is to be noted that the- manner in which this occurs is 
different than in the numerical solution. For the single term approxima- 
tion, the amplitude of the disturbance vanishes, while in the numerical 
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solution, the amplitude of the disturbance stays quite constant as the 
width of the reversed-f low region decreases* The one-harmonic solution is 
displayed on Figs. 3.4a,e to show the differences. Obviously it is 
impossible for the simple harmonic wave to represent the sudden drop into, 
and recovery from, reversed flow which occurs when 3 is near W, consistent 
with our previous remarks about the strengths and weaknesses of the 
approximate method. We will return to this topic when we discuss the 
combined surge and rotating-stall results. 

To summarize, the foregoing discussion of pure rotating stall, as 
well as the results in Appendices A and B, support the following 
conclusions regarding the method of analysis to be applied in the more 
complex combined situation: 

(1) The approximate relation dh/de* = -g seems to be a valid one, 
although no direct proof yet exists. It is also consistent with the 
one-term representation of the solution, and it will be used throughout the 
rest of this report. 

(2) A single harmonic Galerkin method is also a useful approximation. 
It correctly forbids any solution beyond the stability limits of the 
compressor characteristic, and it gives a qualitatively correct picture of 
the nonlinear effect of rotating stall on performance. However, it cannot 
represent the relaxation type of transition which arises in fully developed 
rotating stall, nor can it show effects of diagram steepness. These 
deficiencies are perhaps not serious for the early stages of development of 
a rotating stall transient. In any case, no analytical improvement has 
been found, and we shall adopt the single-harmonic Galerkin method for the 
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remainder of this study. 

(3) The underlying compressor characteristic for this study is the 
axisymmetric one which, as Fig. 3.5 shows, is quite different from that one 
would measure during rotating stall. Accepting the arguments of Ref. [7] , 
we adopt a smooth curve, of which the cubic is a very good model, to 
represent the axisymmetric characteristic. Its only shape parameter is its 
steepness, H/W. 



23 


4. SURGE ANALYSIS 

Governing Equations for Pure Surge 

We turn now to the situation when axial flow coefficient is constant 
around the annulus (not a function of e), but varies with time t. In 
effect, we replace Eq. (3.1) by 

4> ® 0(t) (4.1 ) 

The time-dependence of $ represents surge. As explained and analyzed in 
Refs. [5] and [6], the surge phenomenon depends not only on the compressor 
and its connecting ducts, as rotating stall does, but also on the system 
to which the compressor connects. 

The governing equations (from Ref. [5]) may be expressed as the follow- 
ing pair; the first represents momentum and the second, mass conservation: 

% “ + T(?) - * ($(?)) = o (4.2) 

c d£ c 

* If ■ “T [$(5) “V I - y(?)] = 0 (4 * 3) 

* 4B Z T 

Equation (4.3) assumes a parabolic throttle characteristic, K T ; if the 
throttle can be considered linear, then (bars represent mean values) 

* U [«(5> - T - r- (*<5> - 7>J - 0 (4 * 4) 

C 4B 2 *T 

replaces Eq. (4.3). The dimensionless time and the nondimensional 

effective length of the compressor and its ducting are defined as: 

r = Ut/R ; Z = L /R (4.5) 

^ c c 

Equation (4.3) also implies that the inertia in the throttle duct has been 
neglected. This is generally a very good assumption, because throttle 
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slopes are commonly steep, throttle mass excursions are small, and throttle 
lengths are usually short compared to compressor lengths. 

The nondimensional parameter B, defined in Eq. (1.6), will play an 
important role in the combined problem, as well as in the pure surge case. 
This parameter can be regarded as a measure of the ratio of pressure 
forces to inertial forces for a given rate of mass flow change. High 
values of B thus tend to be associated with system oscillations (surge) 
while low values tend to lead to rotating stall, i.e., operation at a 
stable system equilibrium operating point. 

Without carrying out any detailed calculations the role of the 
different parameters and variables can perhaps be seen most readily if we 
combine Eqs. (4.2) and (4.4) to give the familiar form of a second degree 
ordinary differential equation. This can be regarded as an equation 
describing the behavior of a nonlinear mass/spring/damper system and is: 


d 2 $ 



+ 



4B 2 k_ 


d<) < 

d$" 


d$ 

<3? 


1 


2 2 
U B 
c 


[($ _ $) 



(4.6) 


Several points can be seen from Eq. (4.6). First, it is clear that 
the main part of the "damping" in this system is due to the axisymmetric 
compressor characteristic. The role of the B parameter can be seen to 
affect not only the damping (varying the B parameter can change negative to 
positive damping) but also the stiffness, or restoring term. Finally, even 
though this equation is similar in some respects to that which governs 
rotating stall, one should note that the coefficients are different, and 
the independent (time-like) variable is different, because the basic 
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some respects to that which governs rotating stall, the coefficients are 
different, the independent variable (time) is different, and the basic 
physical situation being described is different* 

We will not belabor the discussion of the basic equations that apply 
to pure surge, since these have been examined in detail in [8] and [9]. 

The two main points that we wish to emphasize as being of importance for 
this study, however, are: 

(1) A compression system consisting of compressor, plenum, and 
throttle is sufficiently general to exhibit surge as well as rotating 
stall. It is therefore presumably adequate for our initial explorations of 
the combined transients, although the ultimate interest is in the 
development of more accurate models of entire engines. 

(2) The axisymmetric characteristic which governs pure rotating stall 
also governs pure surge, in both cases providing the damping effect (which 
can be either positive or negative) of the oscillation. We will find the 
same situation in modelling the combined rotating stall/surge problem. 
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5. DERIVATION OF EQUATIONS FOR GENERAL SYSTEM TRANSIENT 
Pressure Balance of a Compressor 

We turn now to the analysis of oscillations of the compression system 
which are neither pure surge nor pure rotating stall. We expect that such 
motions, which were discussed in the previous two sections, should fall 
out as special cases of the system we will now derive. 

An overall sketch is provided in Fig. 2.1. An irrotational, inviscid 
flow is imagined to proceed from an upstream (atmospheric) reservoir at 
stagnation pressure p T through an entrance duct to the IGV entrance at 0. 
We suppose that <J> , the axial flow coefficient there, can depend on both 
angle 0 around the wheel and time t, although the reservoir pressure is 
constant. As a result of the flow process in the approach duct, if <j> 
varies with 9, then a circumferential velocity coefficient h(£,e) must be 
present at the IGV entrance. It should be noted that the present 
coordinate system is fixed in the laboratory, whereas in Section 3 above 
(as in Refs. [3] and [4]), the system traveled with a supposed permanent 
rotating stall pattern. 

The following averages may be defined: 


around wheel: 


over long time : 



<t>(t,e)de = $(t) 


T 


T __ 

$(t)dt = $ 
o 


(5.1) 


(5.2) 


We recall (Eqs. (2.2, 2.3) that time is measured in radians of wheel 
movement and distance, both circumferentially and axially, in wheel radii. 
Proceeding just as in Refs. [3] and [4], we suppose that any 
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circumferential non-uniformity within the compressor must, by a continuity 
argument and by evidence cited in Ref. [2], go straight through the 
machine. That is, corresponding to Eq. (3.1), we write 

4> = $(5) + g(E,e) ; h - h(?,e) (5.3) 

and note that, by definition, the angle-averaged g must vanish. Also, 
because no circulation can arise in the entrance duct, h must have 
vanishing average as well. Thus, 


2tt 

g(£,e)de = 0 ; 
o 


2it 

h(s,e)de = o 

o 


(5.4) 


For the compressor itself, we recall the general scaled pressure-rise 
formula of Eq. (2.7), and the definition of a in Eq. (3.4). For N stages, 


P F - Pi 1 

— — L - NF - (2$ + 

P D 2 20 5 


(5.5) 


Pressure Balance of the Entrance Duct 

What is wanted is the overall pressure rise, from inlet to exit, 
including the pressure difference associated with the circumferential 
velocity component just ahead of the inlet guide vane. The pressure 
difference from station 0 to station 1 (where the flow is axial) can be 
written as: 


p 0 - p 1 _ 1 „ h 2 

—2 2 K G h 

pU 


(5.6) 


If the IGV entrance is lossless, the entrance recovery coefficient Kq= 1 , 
but if a loss does occur, then Kq< 1 • The exact value of K G probably 


has little effect on the overall results. 
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In the entrance duct upstream of the IGV, we assume irrotational flow, 
so that a velocity potential exists which, though unsteady, will satisfy 
Laplace's equation. Specifically, we define <j> so that 

d> = v/U ; l = u/U (5.7) 

n e 

Accordingly, at the IGV entrance (point 0 on Fig. 2.1), 

4 ) = <H?) + g(£,e) ; U n ) « h(?, 9 ) (5.8) 

no 0 o 


Far upstream, in the reservoir, we take <j> itself to be zero because the 
flow is at rest. 

As a general matter, we must solve Laplace's equation for $ and then 
apply Bernouilli's equation to evaluate pressure at the point 0. For 
unsteady flow, with our particular definitions of variables, the result is 



+ h^ ) + ( <j> ) 

5 o 


(5.9) 


The unsteady contribution, (^)q> is due to unsteadiness in both $ and 
g. The complexity of this term is reduced if we consider a straight inlet 
duct, preceeded by a much shorter potential nozzle, as sketched in Fig. 

5.1. In the duct, of dimensionless length £ j , the angle-averaged 
coefficient $(£) will be constant. Therefore, in view of Eqs. (5.7) and 
(5.8), the velocity potential may be immediately written as 

<t> = (n + 2- I )$(5) + 4>'(5,n) (5.10) 


where $ 1 is a disturbance velocity potential vanishing at n = -£j» and 
giving g and h at point O: 
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(4> * ) = g( 5 ,e) ; (*' ) = h(?,e) 

no 0 o 


The term needed for Eq. (5.11) is simply 


(* r ) = i z If- + (♦') 

5 o 1 d 5 K o 


(5.11) 


(5.12) 


Pressure Balance of the Exit Duct 

In the exit duct, a complicated rotational flow appears even with 
axial OGV's when axial flow varies with 0. We will use an approximation, 
somewhat generalized, developed in [3]. We consider the function P, 
defined below, to be sufficiently small that it satisfies Laplace's 
equation. 


P E 


pU 2 


V P = 0 


(5.13) 


That is, pressure in the exit duct is assumed to differ only slightly from 
static pressure at discharge (p s ). Now, the streamwise Euler equation, 
evaluated at the compressor discharge (point E of Figs. 2.1 or 5.1), gives 
exactly 


( p ) = <♦ r ) = | f + (*'-) 

n e n? o n? o 


(5.14) 


Thus, the potential problem for P is the same as the one already mentioned 
for the minus sign is necessary because, in the exit, n points in 

the streamwise direction, whereas the reverse is true in the entrance. 
Choosing a constant to make P = 0 at the duct exit, rj = fcg, we find 


d$ 

p - V df- 


(5.15) 
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The foregoing simplified analysis assumes that ig is not smaller 
than the distance (conservatively! st which the entrance flow distur— 

bance potential vanishes. If it were much shorter, presumably the second 
term of Eq. (5.15) should be omitted. Therefore, just as in [3], we 
introduce a parameter m in the final result for exit pressure change 


P S~ P E 

P u 2 


d$ 


= < p >* = - Cm-1 > C ♦') 


E d? 


S O 


( 5 . 16 ) 


»here m = 2 would refer to a "long enough" exit duct, and m = 1 would 
refer to a very short one. 

Net Pressure Rise to End of Compressor Exit Duct 

Between the upstream reservoir (p^) and the discharge from the exit 
duct (pg) is the zone where circumferential pressure variations may 
arise. Further downstream are the plenum and throttle, which we assume to 
contain only axisymmetric disturbances. We can now form the net pressure 
rise in the zone of possible angular variations by combining Eqs. (5.5), 
(5.6), (5.9) and (5.16). After introducing Eqs. (5.3) and (5.12), the 

result may be written 

Po-Pn 


” SJ r m w-i* 2 ’ - (»! + 7 + V f 


pU 


-h (2 *?n **on’o ' 2 " - K G ,h 


(5.17) 


Following [3] we identify the left side, which is the total-to-static 
pressure rise coefficient, as <?(?) (which also appears in Sections 3 and 4 
above), and the first parenthesis on the right as the quasisteady, axisym- 
metric compressor characteristic, also familiar from Section 3; we notice 
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that for steady flow with no disturbance, all other terms on the right 
vanish, leaving the obvious identity ¥ = il> c . 




pU 


H c (4>) 2 NF(*) - - 


(5.18) 

(5.19) 


The second parenthesis on the right of Eq. (5.17) is the effective 
flow-passage length through the compressor and its ducts, measured in radii 
of the compressor wheel. This has been denoted as t c and thus 



As explained in [3], if the compressor lag t is considered to be purely 
inertial, one could write 


_ 1 _ 

a 


(2NL R > (i) -V 

cos y 


(5.21) 


The first factor in parentheses is simply the axial length of the 
compressor, being the axial length of a row; the factor 1/R puts the 
distance into wheel radians; k is a factor which could account for 
interrow spacing; and cos^y accounts for the tortuous flow path through 
the compressor, with y being the effective stagger angle of the blades. 
Thus, the effective (inertial) path length for the system is longer than 
the axial length. The total path length % c has already been mentioned 
in Section 4 above. 

The last term of Eq. (5.17) will be neglected, in effect assuming 
that the recovery coefficient K G at the IGV entrance is 1. Finally, 
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therefore, we may write Eq. (5.17) in the form 

,,t ’ ( *;v -*«t 

In this we have used Eqs. (5.3) and (5.11) to express the argument of 
in terms of Z ' : 


= $ + ( $ ' ) 

n o 


(5.22) 
♦ c 


(5.23) 


The next step will be to place Eq. (5.22) in the context of the complete 
compression system, with plenum and throttle, as sketched in Pig. 5.1. 


Before doing so, we simplify Eq. (5.22) by introducing a generalization of 
an approximation familiar from Section 3. 

An Approximation of the dh/de = -g Type 

The disturbance potential * satisfies Laplace 1 s equation 


t>. + <(> ' =0 

00 nn 


(5.24) 


Its value and normal derivative <(> 1 ^ at the compressor entrance are 
needed, which implies that Eq. (5.24) must be solved as part of our 
analysis. We would like to avoid the difficulties of such a procedure, if 
possible. 

We know that g is a periodic function of 0, and, by definition, must 
have a vanishing average over 2ir . Therefore, it has a Fourier series and 
the solution of Eq. (5.24) is 


CO 1 

I — e n (a sin n0 + b cos n0 ) 
n =1 n n n 


T) < 0 


(5.25) 


vanishing at r) = 


- 00 . 
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Now, if only the first term of Eq. (5.25) were kept, one would 
conclude that 



- -(<!>' ) 
00 o 


(5.26) 


which, in view of Eq. (5.11), is precisely the familiar dh/d0 = -g 
relation. Thus, it will be expedient to use the value of <j> 1 at t] = 0 as 
dependent variable, using the notation Y for simplicity: 


(r) = Y(K,9) ; (♦') = -Y 

0 T| o 00 

Now Eq. (5.22) takes the form 

d$ 1 

^ v c 99 c d£ l 2a £99 999 


(5.27) 


(5.28) 


Eq. (5.27) implies 

h . Y 8 ; g - -Y e# (5.29) 


Therefore, if Y is simply required to be periodic in 0 , as it obviously 
must be to describe a physical quantity in the compressor, 

Y( £ , 0+2tt ) = Y(£,e) (5.30) 

then the cyclic integrals of g and h vanish, 


2tt 

g( 5 , 0 )de = 0 

0 


2tt 

h( 5 ,e)de = o 
o 


(5.31) 


as we know they must do; in the case of g by definition, and in the case 
of h because no circulation can arise in the inlet as we have described 
the flow there. It is also obvious from the general solution of Eq. 
(5.25) that the cyclic integral of <j> 1 itself must vanish, so we know that 
Y should have the property 
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2ir 

Y(?, 9 )de = 0 (5.32) 

0 


The idea described in Appendix A to improve the accuracy of the 
dh/d9 = -g approximation for rotating stall can also be adapted for the 
present more general case. Corresponding to Eq. (A.6), we write 



A 


-Y 


00 


■ A A A A 

( <|> „ ) = AY - BY 

0 0 6 9 9 9 


(5.33) 


Just as in Appendix A, consistency between Eqs. (5.25) and (5.33) can be 
arranged for two harmonics. Choosing the first two, one finds A = 2/3 

A 


and B = 1/3, as in Eq. (A. 8). Corresponding to Eq. (A. 9), one would 
find the following equation when Eq. (5.33) is used to modify Eq. (5.22) 


, * * d$ 1 ^ 1 ^ 

¥ = ($ - Y ) - l -tt - mAY +(- L +mB)Y + —— Y 

c 00 c £ a £00 2a 900 


(5.34) 


As described in Appendix A, the change is a matter only of coefficients. 
For present purposes, we will omit this refinement, but keep it in mind 
for future studies. 

Overall Pressure Balance of the Compression System 

We now return to Fig. 5.1 and account for pressure changes downstream 
of the compressor exit, in the plenum and in the throttle. We will derive 
relations somewhat more general than those already mentioned in Section 4 
above . 

The plenum will obviously eliminate any spatial variations; thus, it 
will receive mass at a rate pUA c $(?;), and discharge mass at a presumably 
different rate pUA c $ T (H;). Duct areas at entrance and discharge of the 
plenum would both presumably be different from the compressor area A c ; the 
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difference is accounted for in the definition of as in Ref. [5]. 

If $ and are different, then mass must accumulate in the plenum, 
changing the density there. As discussed in [5], we can take this change 
to be isentropic for the situations under consideration. 

We assume that the compressor exit discharges as a free jet into the 
plenum, so that its pressure is p s (£). The isentropic relation between 
density and pressure then provides that the rate of density change in the 
plenum is the ratio of dp s /dt to the square of sound speed, a s . The 
rate of increase of stored mass in the plenum volume (Vp) is thus propor- 
tional to the dimensionless time derivative of the pressure-rise 
coefficient ¥(£), defined in Eq. (5.18). The balance of entering, leaving, 
and stored mass of the plenum can then be written in the following 
convenient form (see Eq. (4.3)): 

- i " 1 ' 1 - V ' 11 (5 - 35) 

> 4B 


The plenum discharge coefficient is of course also the throttle 
flow coefficient. It must be related to the pressure difference applied 
to the throttle, and the throttle performance characteristic, in order to 
complete our overall pressure balance. Generally, these relations would 
be selected to properly represent an engine application. Here, we 
consider that the plenum discharges without loss into a throttle duct of 
length L T> in which there is a throttle with a pressure-flow coefficient 
characteristic given by 


= F ( $ ) 
2 T ' T ' 
pU 


(5.36) 



36 


The throttle discharges to a reservoir (atmosphere) at the same pressure 
prp as that which earlier supplied the entrance to the compressor (see 
Fig. 2.2). Therefore, the pressure difference Ps~Pt balances the 
throttle loss and also provides any acceleration of the mass in the 
throttle duct. In coefficient fora, 


d* 

»(?> - vv + Sir 


(5.37) 


In the present paper, we will take the throttle duct as short enough 
to justify neglecting Jl T . This is generally quite a good assumption for 
realistic throttle ducts. As discussed in Section 3, two typical 
possibilities for the throttle characteristic, F T ($ T ) are: 


parabolic: Y = 

linear: Y = Y + k T ($ T ~$) 


(5.38) 

(5.39) 


In any case, having neglected d T , we write that $ T is simply the 
inversion of F T : 

^T = Fj 1 (^ ) (5.40) 


In principle, Eqs. (5.35) and (5.40) complete our system of equations 
for disturbances of the compression system. We imagine that $ T may be 
eliminated between these two equations, and the resulting equation 
connecting ¥ and $ can be combined in some way with Eq. (5.28). If so, 
there are only two equations, so far, for the three unknowns ¥(£), $(£), 
and Y( £ , 0 ) • The third equation comes from realizing that Eq. (5.28) 
involves terms which are functions only of time. If that equation is 
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integrated over a cycle of 0, the following relation results: 


/x d$ 1 

y( 5) + a, — = — 

c 2tt 


o'*o<‘ - X 0 8 ,de 


(5.41) 


This integral relation is, in effect, the needed third equation. 

Equations of a General Disturbance 

For convenience, we assemble the three pressure-balance equations we 
will use in subsequent analysis, giving them new numbers. 

(5.42) 

yy ^ za t,yy yyy 

’2tt 


+ *C df - 'c 1 ' - V - " Y 5 + k (2Y ;e« + \ee> 


’ ( 5> + t c«-ir 


, (« - y >de 

o C 00 


d¥ 1 -1 

^cd? = -2 [$(?) - F T 
s 4B 


(5.43) 

(5.44) 


where partial differentiation is indicated by subscripts £ and 0. 

Pure Rotating Stall and Pure Surge as Special Cases 

It is instructive to see how Eqs. (5.42-5.44) specialize for the pure 
rotating-stall and surge motions outlined in Sections 3 and 4 above. 

Surge is the simplest; in that case there is no 0 variation, and Y(£,0) 
may be set to zero. Eqs. (5.42) and (5.43) are then the same equation, 
identical to Eq. (4.2). Eq. (5.44) is, of course, identical to Eq. (4.3) 
for a parabolic throttle line. 

For pure rotating stall, the plenum pressure, and the angle-averaged 
through flow are constant with time; therefore, dy/d£ and d$/d£ both are 
zero. The right side of Eq. (5.44) vanishes, becoming equivalent to a 
statement that $ = $»p, and that the operating point must lie on the 
throttle characteristic (e.g., Eq. (5.38)). Y can be represented as a wave 
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traveling at speed f; that is, a function of a single variable: 

Y(5,0) i Y(0*) ; 9* = 0 - f£ (5.45) 


The variable 0 1 is the angle measured in the disturbance frame, and was 
defined in Section 3. Eqs. (5.42) and (5.43) become, after recalling the 
definition of X _rom Eq. (2.3): 


, d 3 Y c dY r 

X + mf - [Y - t|i ($ 


d0’ 


d0 


d 2 Y 

d0* ; 


)] = 0 


Y = 


2tt 


2tt 


2 

Y _ )d0* 
0 C d0* 2 


(5.46) 


(5.47) 


Next, we recall Eq. (5.29) and see that now 




-g 


(5.48) 


This last equation converts Eq. (5.46) precisely to the governing equation 
of pure rotating stall, Eq. (3.2). Eq. (5.47) is converted to 
1 


Y = 


2ir 


2-j: 

i|> ($ + g)d0 ; 
0 c 


(5.49) 


which the definitions of 6 and G in Section 3 show to be identical to Eq. 
(3.12) for performance effect of rotating stall. Thus, the reduction of 
our general disturbance equations to pure rotating stall is demonstrated. 
Application of One-Term Galerkin Procedure 

To treat a general disturbance which may have both angular variation 
like rotating stall and time-dependent mean flow like surge will require 
solution of the complete system of Eqs. (5.42-5.44), which is third order 
in angle, but only first order in time. This circumference suggests that 
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a Galerkin treatment of angle dependence would be helpful to eliminate the 
angle dependence and reduce the system to three ordinary differential 
equations, with time as the only independent variable. 

Section 3 above indicates how to proceed with the Galerkin procedure. 

We have two issues to settle before beginning, however. First, should we 
imitate Section 3 closely by working in a traveling frame (0*), or should 
we stay in the "laboratory frame" (0)? Study shows there is no advantage 
and many conceptual difficulties in the former system, while the latter is 
perfectly general, depending only on the obvious consideration that 
physical quantities must be periodic in wheel angle whether or not a perma- 
nent wave is in progress; therefore, we will apply the Galerkin procedure 
in laboratory coordinates. 

Second, shall we be satisfied with the one-term harmonic representa- 
tion of angular variation? For the present study, there seems no good 
alternative if we want to keep a relatively simple form. Even if we found 
a better functional form to represent the rotating-stall profiles of Figs. 
3.4a-f, it is not clear that that would be helpful for general disturbances, 
some of which might be weak transients which actually would be well 
represented by a simple harmonic wave. Although we might consider 
extending the analysis to higher harmonics, experience with a two-harmonic 
approach (for pure rotating stall) suggests that such an extension would 
have to be carried to high harmonic order to be valid. For the present, we 
will accept the accuracy limitations inherent in a one-harmonic model. 

Recalling Eq. (3.17), we represent Y as follows: 


Y = WA(£) sin( 0 - r(?)) 


(5.50) 
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We include phase angle r(£) in this equation, because we must not expect 
nodes of the disturbance always to remain in the same angular location. 


We substitute Eq. (5.50) into Eq. (5.42), and find that the residual 
(analogous to Eq. (3.18)), moments of which will be made to vanish, 

W L " _ w I ’ r + * c f§'“ + WA sin(e-r))] - (m+1) || sin(e-r) 

. r/ , JL\ dr 1 , „ 

+ [<nH 7 ) dF ■ 2l lA cos(0 - r) 


is 


(5.51) 


where Y , $, A, and r are to be functions of 5 alone. 

We now set the moments of Eq. (5.51) equal to 0. The integral of L 
must vanish; this accounts directly for one of the basic equations, Eq. 


(5.42) : 


_ 1 _ 

2ir 



+ WA sins = y + $, 


d$ 
c d£ 


(5.52) 


The sin(Q-r) moment is 


1_ 

irW 


2tt . . . 

sins i|> ($ + WA sins)ds = (m + -) ^ 
0 c a d£ 


(5.53) 


and the cos(e-r) moment is 


1 

irW 


f 2ir 4 j . 

^ coss 't' c (^ + WA sins)ds = -t(m + — ) 


(5.54) 


Without yet specifying the characteristic function i|j cf we see that 
if it is a regular function of $ , the integral on the left must vanish 
when carried out over a complete cycle. Thus, very generally, it is clear 
that dr/dj; must be a constant if A is not to be zero. Calling that 
o» 


constant f, 
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r = f • l 
o 


f = -1Z1 — 

o ~ 1 + ma 


(5.55) 


Actually, f Q is the propagation speed found in [3] for small-disturbance 
rotating stall (which is a simple harmonic wave). We conclude that, even 
though the amplitude A varies with time within the context of the 
one-term approximation , nodes of the wave of angle dependence must travel 
at a constant speed comparable to that of pure rotating stall. 

We are left with Eqs. (5.44), (5.52), and (5.53) to determine the 
unknown functions of time , Y , $ , and A. 

Final Simplified Equations 

It remains to specify the axisymmetric characteristic; we adopt the 
cubic function defined in Eq. (3.21), making the substitution 


•^7 - 1 = 7T~ 1 + A sin£ 


W 


W 


(5.56) 


It is obvious that if can be expressed as a sum of powers of <j> , only 
even powers of A sin^ will contribute to the integral in Eq. (5.52), and 
only odd powers in Eq. (5.53). Inspection of those equations then shows 
that A will enter as the square, no matter what the particular form of 
\|> c . For convenience, we define a new variable J to replace A: 

J(?) = A 2 (?) (5.57) 

and we note that J must always be positive. Carrying out the indicated 
integrals, we find the governing equations to be 
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dJ 

at 


J[i 




3a H 
( 1 +ma)W 


(5.58) 

(5.59) 

(5.60) 


These are our final equations for instantaneous values of flow coefficient 
($), total-to-static pressure rise (¥), and squared amplitude of angular 
variation (J), all as functions of time ( 5 ). 


1 
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6. DISCUSSION OF SIMPLIFIED EQUATIONS 
Pure Modes and Their Growth 

Our intention is to use the foregoing Eqs. (5.58-5.60) to compute the 
consequences of an initial disturbance of a compression system. Before 
doing that, we should examine these equations to discover what we can 
about their general features. 

First, we observe how these equations describe pure surge or rotating 
stall. In the case of pure surge, J = 0, and Eq. (5.59) becomes simply 
Eq. (4.2) with cubic \[> c (Eq. (3.13)). Eq. (5.58) is the familiar Eq. 

(4.3), in effect. In pure rotating stall, amplitude J is constant, and 
the bracket on the right of Eq. (5.60) must vanish, establishing an 
,, equilibrium ,, J e : 

J = 4[1 - (J - I) 2 ] (6.1 ) 

e w 

In this case, $ will also be constant because there is no surge, and we can 
identify ($/W) - 1 as the quantity B introduced in Section 3 to describe 
average flow coefficient, and show that J e is simply the A*^ for pure 
rotating stall in Eq. (3.20). We may also observe that when this value of 
J is inserted into Eq. (5.59), and d$/d£ is set equal to zero, that 
equation reduces to the one-harmonic Galerkin solution for performance 
effect (6) in rotating stall, presented in Eq. (3.22). 

Our equations appear to permit the existence of pure modes; that is, 
either surge or rotating stall without the other. However, we should 
question whether such modes could evolve from initially infinitesmal 
disturbances. Clearly, such a process is possible for pure surge; if J is 
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zero it will remain zero, according to Eq. (5.60), and the other two 
equations will deal with the evolution of pure surge disturbances, as 
described in Refs. [5] and [6], 

It appears to be impossible, however, for rotating stall to evolve 
without producing at least some disturbance of $ or if J changes, then 
Eq. (5.59) says that K or } must also change. A disturbance of $ induced 
in this way might of course be a non-cyclic transient, not to be 
identified as "surge". Small values of B are associated with rotating 
stall in preference to surge (Ref. [5]). The present Eqs. (5.58), (5.59) 
suggest that if B is small enough, and the throttle is steep enough, $ 
could remain essentially constant while J grows. The solution of Eq. 


(5.63) in 

that case 

is 


J 


1 


J 

e 

i 

d)| 0 

+ 

T— 

1 ) exp 

' 3aJ H/W 1 
e 


4 ( 1 +ma ) ^ 


where J 0 is an assumed initial disturbance. The essence of the solution 
is that J grows steadily from J Q to its final, fully-developed pure 
rota ting-stall value of J e . 

Figure 6.1 shows this behavior for the compressor geometry discussed 
in Section 4. For reference the relevant parameters are: m = 1.75, H = 

0.18, W = 0.25, 1/a = 3.5. The calculations have been done for three 
different values of J e /J Q : 10 + 2 > 10 + ^, and 10 + ^. If we consider 

the transient to occur at $ = 0.25 (i.e., 8=0), then this implies initial 
rotating stall amplitudes of 0.2, 0.02, and 0.002, respectively. 

The graph shows the normalized rotating stall amplitude A/A e (= /J/J e 
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versus time, with time being plotted in rotor revolutions (1 rev = 2 ir£). 
It is clear from the figure that the growth is dependent on the initial 
conditions. This can be seen from Eq. (6.2) which, if J Q «J e , for 
small times can be written as 
3a J H/W 

J - J = J (- 7 — 77 - r)E + ... 

o o 4a(1+ma) s 

i.e., the increment of J is proportional to the initial value. Once the 
value of J has grown so that the nonlinear aspects of the process become 
important, however, the growth rate to the final limit cycle amplitude 
becomes independent of the initial conditions. This can be seen by 
noting that if the curves for J e /J c = and were translated 

(in time) such that they matched with the curve for J = 10 4 at the point 
where J/J e = 0.2, say, the three curves would be virtually identical 
all later times. This type of behavior, which is typical of nonlinear 
systems, will also be encountered later, when we discuss the results for 
the combined rotating stall/surge transients. 

Eq. (6.2) says that, for any given value of £$, the time needed to 
approach J within a closeness defined by 


J -J 
e 


J 

e 


e 


(6.3) 


is 


K 


e 


1 +ma 

3a( 1-6 2 /w 2 )h/w 



(6.4) 


Thus, if J Q «J e , the time of final approach is proportional to in J Q 
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Also, the time for rotating-stall development is short for steep diagrams 
(large H/W), and for throttle settings just inside the stall zone (6 less 
than but close to 1 ) . The fastest approach would occur if the compressor 
were already (transiently) in the middle of the steady-state rotating stall 
zone (8 = 0) before rotating stall first appeared. 

Nature of the Coupling Process 

We have seen that Eq. (5.60) describes how angular— disturbance 
amplitude J grows at a rate which initially depends only on the 
disturbance itself, but then, in the manner of a chemical reaction, tries 
to approach an equilibrium described by Eq. (6.1). In general, the flow 
coefficient will change with time, so that during a coupled oscillation, 
J e U) will be a moving target, so to speak, for J. Whether J actually 
achieves and stays on that target, that is, quasisteady equilibrium with 
$(?), will be known only from simultaneous solution of the equations, 
giving in effect a trajectory in the three-dimensional space , $, J. 

Inspection of Eq. (5.60) also shows that excursions of $ away from 
the central point $ = W will diminish the rate of change of J. Thus, if 
we suppose a throttle setting of $ =W (i.e. , 6 = 0), then the presence of 
surge-like variations of $ will tend to suppress circumferential 
variations. Thus, if B is small, Eq. (5.58) implies that swings of $ will 
tend to be small, which, relatively, would encourage circumferential 
variation, according to Eq. (5.60). 

Turning to Eq. (5.59), we see that the presence of a circumferential 
angular variation in velocity will diminish the term in which J appears 
(recalling that J must always be positive). That term describes the 
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central slope of the characteristic, at the midpoint of Fig. 3.2, thus 
reducing the basic amplifying effect of the characteristic in the stall 
region. We can conclude that this variation (i.e., the presence of a 
rotating stall-like disturbance) would tend to inhibit surge-like 
disturbances of $. 

The foregoing arguments, which are certainly in qualitative agreement 
with the many experimental observations that have been made of these 
phenomena, only suggest how "surge-like" and "rotating stall-like" 
disturbances are coupled through Eqs. (5.59), (5.60). The coupling is 
mediated by Eq. (5.58), and one must carry out trajectory calculations to 
see the actual consequences of a given initial disturbance. 

Small Surge Disturbance of Rotating Stall 

We know that our equations permit both pure surge and pure rotating 
stall. We might go further and ask whether those "pure" motions are 
stable to disturbances of the other family. First, we imagine a case of 
fully-developed rotating stall subject to weak surge-like disturbance. 
Assuming an infinitely steep throttle line (kj = <*> in Eq. (4.4)), and 
eliminating Y between Eqs. (5.58, 5.59) by differentiation, while holding 
J constant, one finds a harmonic-oscillator equation for $(£), with the 
damping term 


damping : 



2 „ „ 

b _L iii 

2 J W dr 
W s 


(6.5) 


Damping would be negative, if J were zero. However, if J is fully developed, 
then J is actually 4(1 - B 2 /W 2 ), which makes Eq. (6.5) positive: 
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damping : 


2 r i B_, 1 d$ 

2 11 ~ w 2 J W d C 


and we conclude that small-disturbance surge-type oscillations tend to have 
positive damping in the presence of "equilibrium" rotating stall. There- 
fore, any surge-like disturbance of pure rotating stall would therefore 
decay, which might be inferred from the fact that the relevant compressor 
curve is now the negatively— sloped fully developed rotating-stall curve. A 
careful stability analysis of Eqs. (5.58-5.60) should be carried out; quite 
possibly the matter is more complex than the foregoing argument implies. 
Small Angular Disturbance of Surge 

Next, we suppose that an equilibrium surge oscillation is in 
progress, and ask what happens when a small rotating-stall-like 
disturbance (J Q ) is introduced. Eq. (5.60) will govern the growth or 
decay of J. We note that $ will vary during the surge oscillation so that 
the bracket of Eq. (5.60) may be either positive or negative, and it is 
therefore not obvious whether, on balance, J will grow or not. 

To estimate the effect of $, we recall the pure surge equation, Eq. 
(4.6), and set k T = a> (steep throttle line). Then, the surge equation 
takes the same form as Eq. (3.9) for rotating stall, with the time 
variable 5/(2Bi c ) playing the role of 0*. We can therefore adapt the 
Galerkin solution for rotating stall (Eq. (3.21)) to represent surge: 


$ _ $ 



sin 



( 6 . 6 ) 


Identifying ($/W) - 1 as 6/W and introducing Eq. (6.6) into Eq. (5.60) 
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gives, after simplification, 


1 dJ = 3aH/W . el 

J 1+ma 1-1 2' 

w 



- sin 

VT 



+ 2(1 cos2 ffiT " 7 11 

W c 


(6.7) 


If we then neglect J within the bracket on the ground of small 
disturbances, the foregoing equation can be integrated to give the level 
of J predicted to exist at any time £ during the disturbance: 


J 


exp 


3aH/W 

1+ma 


2B£ 
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[-(1 



2BSL + 
c 





2B l 

c 


+ 


(1 - sin2 

W 



( 6 . 8 ) 


where a constant of integration is included. The question is, 
whether J is greater or less than its initial value J Q . 

If the curly bracket of Eq. (6.8) is positive, J>J 0 , while, if it 
is negative, J<J Q . Examining that expression in detail, one sees that 
the first term in the square bracket, being negative and increasing, 
represents a tendency for J to be progressively smaller than J Q for 
large time. The oscillatory harmonic terms could reverse that tendency 
for early or moderate times. 

Fig. 6.2 illustrates the possibilities for the special case of B = 0. 
The square bracket, the sign of which indicates whether J is greater or 
less than its initial value J c , takes the form 
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* s a r 

growth factor: [- ■ + sin2 ] (6.9) 

c c 

The first and second terms of this expression are plotted separately on 
Fig. 6.2. A series of lines are shown to represent the first term, because 
is arbitrary. Where any one of these lines cross the sine wave is a 
possible starting time for the disturbance. What happens after the start- 
ing time clearly depends on the phase of the sine wave (surge) when the 
angular variation starts. If the start is at £ = 0, then the shaded area 
above line (1) shows the duration and intensity of the excursion of J. The 
duration is less than a quarter period of the surge oscillation. Whether 
the amplification rate is sufficient during this period for rotating stall 
to develop cannot be found from this small-perturbation solution. 

If the intersection point is moved along the sine wave, the time 
available for amplification shortens, until line (2) is reached for which 
there can be no growth at all. Not until line (3) is reached does 
amplification again become positive. The greatest tendency for rotating 
stall to grow is found for line (4). 

In order better to visualize these results, we sketch on Fig. 6.3 the 
cubic compressor characteristic, and overlay a surge cycle sketched as a 
counterclockwise circle centered at B = 0. We indicate by dark lines on 
the circle certain phase zones. If an angular disturbance initiates in 
one of those zones it will, at least briefly, amplify. Outside those 
zones, amplification of a weak angular disturbance cannot occur. Points 
are labeled on Fig. 6.3 to show correspondence with the lines of Fig. 5.1. 
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The foregoing discussion indicates that it is very important to know 
the precise character of an initial disturbance, and the timing of its 
introduction, whether it is velocity-like or pressure-like, and its 
magnitude. All those factors could have a profound effect on the ultimate 
trajectory of the disturbance. 

Summary of Qualitative System Behavior 

Before showing the numerical results for the trajectories in $ , y , A 
space, we can summarize the conclusions that we have drawn so far about 
the system behavior: 

( 1 ) Pure modes of either surge or rotating stall are permitted as 
permanent oscillations. 

(2) Pure surge can evolve from an initial weak axisymmetric disturbance. 

(3) In principle, pure rotating stall cannot evolve directly from a weak 
angular disturbance; some axisymmetric disturbance will be induced. 
Small B-parameter and steep throttle characteristic would favor growth 
into pure rotating stall. 

(4) The growth rate of an angular disturbance (towards a limit cycle 
value) depends initially on the strength of the initial disturbance, 
but becomes independent of this when the nonlinear aspects of the 
process become important. 

(5) Angular disturbances will grow toward an equilibrium value depending 
on the instantaneous flow coefficient; if the latter is also changing 
with time, angular-disturbance amplitude "chases a moving target". 

(6) Excursions of flow coefficient will tend to suppress angular 
variations, and the presence of angular variations will tend to 
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suppress surge. 

(7) Any small surge-like disturbance in the presence of fully-developed 
rotating stall will tend to decay. 

(8) In the presence of fully-developed surge, a small angular disturbance 
may decay or may grow (but only during a fraction of one surge cycle), 
depending on the phase of the surge cycle in which the disturbance 
originates. 
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7. NUMERICAL RESULTS FOR GENERAL POST-STALL TRANSIENTS 

The set of coupled ordinary differential equations given as Eqs. 

(5.58), (5.59), and (5.60) have been solved numerically for a 
representative set of parameters to show some of the basic features of the 
phenomena. Specifically, calculations have been carried out to illustrate 
the effects of the non-dimensional parameters B and i c as well as of the 
initial conditions. In view of the approximate nature of the approach 
taken, it is to be emphasized, however, that no extensive parametric study 
has been performed. 

In the calculations, the axisymmetric compressor characteristic has 
been taken to be the same cubic curve that was described in Section 3, and 
this is shown for reference as the solid line in Fig. 7.1. Except as noted 
otherwise, the throttle line used is the parabola that is shown as the 
dash-dot line in this figure. As mentioned in Section 3, these curves are 
representative of a three-stage, low-speed compressor, such as that used in 
Ref. [5]. 

Unless specified, the initial values of non-dimensional pressure rise 
and flow coefficient are taken to be the values at the peak of the curve, 
i.e., $ = 0.5 and ¥ = 0.66. At that condition, we shall imagine that a 
certain level of circumferential nonuniformity defined by A(0) appears at 
time 5=0*. The system of equations (5.58-5.60) governs the ensuing 
transient which, owing to coupling, may show either temporal (surge-like) 
or circumferential ( rotating-s tall-like ) variations. 

*Recall that the nondimensional amplitude of the 6-dependent velocity 
nonuniformity is AW, and A = /J. 
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Effect of B Parameter 

We start by demonstrating how our equations give a more or less 
familiar effect, that of the B parameter. (This has previously been well 
documented [5,6] as far as the overall system response is concerned.) 

Figure 7.2 shows the transient system response in terms of the overall 
pressure rise (atmosphere to plenum) and annulus averaged flow coefficient 
for three different values of B; 1.0, 1.4, and 2.0. Also indicated in Fig. 
7.2 is the axisymmetric curve (the solid line) and the calculated rotating 
stall curve, which is shown dashed. In these calculations, the value of 
8.0 and the initial rotating stall amplitude A(0) was one percent 
of the average velocity, i.e., the initial non-dimensional amplitude of 
velocity perturbation was 0.005. 

It is seen that the higher values of B cause greater excursions in 
axial velocity. At the highest value, this would lead to surge cycles. 

For the lower values, the eventual result of the initial transient is 
operation at a new equilibrium point on the rotating stall curve. 

As mentioned, the picture shown in Fig. 7.2 is a somewhat familiar 
one, but is nevertheless remarkable in that it follows from the present 
general transient equations. Figure 7.3 presents a new, more 
detailed view of these same transients. This latter figure shows the 
rotating stall amplitude, A(£)W (= /JW) , versus annulus-averaged flow 
coefficient ($) for the same three values of B, and the same value of £ c 
and initial conditions. Also indicated in the figure is the locus of 
steady-state rotating-stall amplitudes for different values of $; in other 
words, the amplitude one would have in steady-state operation at that 
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particular mass flow. This is just the value of amplitude that would be 
predicted from Eq. (3.21). For the cubic characteristic used, the steady- 
state amplitude has a maximum at $ = 0.25 and goes to zero (i.e., 
steady-state rotating stall ceases) at $ = 0.5 and $ = 0.0. This curve of 
amplitude gives the "equilibrium" toward which rotating-stall-like 
disturbances must tend. 

An obvious feature of these transients is that over the first part, 
for $ greater than roughly 0.25, the amplitude during the transient is 
substantially less than that during steady-state operation. However, for 
lower flows than this, the amplitude can actually exceed the steady-state 
value. In addition, as seen in the curve for B = 2.0, the rotating stall 
persists into the reverse-flow regime during the transient. This can be 
contrasted with steady-state operation in which, for this compressor, 
rotating stall cannot exist when flow is reversed. 

It should be noted that when the transient is of the surge type, 
calculations have only been carried out for part of the initial cycle. The 
reason for this is that the stall-cell amplitude decays to a very low 
value during the extended operation at reverse flow. Thus the "initial 
condition" for the part of the transient when the flow through the compres- 
sor accelerates from zero back to unstalled operation is essentially that 
of zero stall-cell amplitude. Under these conditions, as discussed in 
Section 6, there will be no subsequent growth of rotating stall, and the 
transient will proceed as pure surge. 

It does not seem physically realistic, however, that there would be a 
lower level of flow asymmetry during a surge cycle than during operation 
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P^ior to stall, and thus this part of the calculation does not appear to 
be a valid description. It might be more correct to say that there is 
always some minimum disturbance amplitude present, able to provide an 
"initial disturbance" at any time. Although there is certainly no bar to 
including this in the numerical integration routine, in view of the lack 
of information about this point, we have not elected to do this, but 
rather to restrict computations to the first part of each transient. We 
will return to a discussion of this point in the next section. 

The basic feature shown in these calculations is the growth and decay 
of the rotating stall, either to zero or to the appropriate "stagnation" 
value, as the mass flow varies. That this process must occur has been 
known qualitatively for some time. As far as the authors know, however, 
this is the first time that it has been calculated, in even an approximate 
manner , using the equations of motion for asymmetric flow in a compressor. 
Effect of Initial Conditions 

In addition to the B parameter, the effect of the initial level of 
disturbance A(0)W was also examined. Results are shown in Fig. 7.4, which 
shows the locus of the compressor operating point in the AW, $ (rotating 
stall amplitude/flow coefficient) plane. The calculations shown are for B 
= 1.0, l c = 8.0, and for three different initial values of rotating-stall 
amplitude: 0.05, 0.005, 0.0005. These correspond respectively to initial 

velocity perturbations of ten percent, one percent, and one— tenth percent 
of the mean axial velocity. It can be noted that the last of these is 
roughly the size that would be expected in a good wind tunnel test section 
and the first appears to be larger than the long-wavelength velocity 
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disturbances that are seen prior to rotating stall. Thus, our range of 
initial disturbance amplitude conditions encompasses essentially all 
physically plausible situations. The value that we regard as our "base 
case", the one percent level, might be a reasonable estimate for the 
magnitude of velocity non-uniformities to be expected prior to rotating 
stall, but there seems to be little data on this point. 

In Fig. 7.4, the locus of steady-state rotating stall amplitude is 
again indicated by the crosses. The eventual equilibrium point in 

rotating stall, at $ - 0.35, ¥ = 0.455, is also indicated as point E. 

The effect of the initial condition is strongly evident early in the 

transient, since the growth rate is proportional to the amplitude of the 
perturbation. After the stall amplitude has reached a certain size, 
however, on the order of 0. 1-0.2, the nonlinear effects become important, 
and the slope of the trajectory does not depend strongly on the initial 
conditions. Although there are some differences in the extent of the 
excursion in axial velocity between the three cases, the behavior once the 
trajectory approaches the vicinity of the steady-state locus is fairly 
similar for the three cases. As discussed in Section 6, this relative 
independence from the initial conditions during the approach to the 
equilibrium state is typical of nonlinear systems of the type which we are 
examining. 

The effect of the initial condition can also be shown in a $ , ¥ 
representation, as in Fig. 7.5. Here, only the curves for initial 
amplitudes 0.05 and 0.0005 are shown. It can be seen that even though 
there are two orders of magnitude of difference between the initial 
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amplitudes, the trajectories of the plenum pressure versus flow 
coefficient are not at all dissimilar; they both show roughly the same 
excursion in mass flow before ending up at the eventual equilibrium point 
in rotating stall, denoted on the graph by E. 

We can also look more directly at the influence of initial conditions 
on the compressor pumping characteristic, i.e., the compressor output. 

This is found by subtracting the momentum changes in the compressor duct 
from the overall (atmosphere to plenum) pressure rise. The compressor 
pumping characteristic is given by 

Nondimensional pressure rise ® \|j ($) + 6 (7.1) 

where now 6 is a function of time. The result of doing this is seen in 
Fig. 7.6. This figure presents the instantaneous compressor pressure rise 
versus the instantaneous annulus-averaged axial velocity parameter, for 
the same conditions as in Fig. 7.5, namely B = 1.0, = 8.0, and 

initial amplitudes A(0) = 0.05 and 0.0005. 

The trajectory for 0.0005 can be seen to follow the axisymmetric 
characteristic until near $ ~ 0.25 (the point of maximum steady-state 
rotating stall amplitude), and then depart from that curve and thereafter 
follow the steady-state rotating-stall curve quite closely. The 
trajectory for 0.05 departs from the axisymmetric curve virtually from its 
initial point (showing the influence of nonlinearity) and tracks 
essentially along the steady-state rotating-stall curve once it has passed 
through for the first time. The eventual equilibrium point is again 
indicated in the figure as E. 
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Calculations have also been carried out for other values of B to 
examine the effect of initial conditions. These confirm, in general, the 
trends shown in the preceeding figures, namely that over a fairly wide 
range of initial conditions, the results do not depend strongly on the 
precise values used. There are, however, situations in which the initial 
condition can make a great deal of difference in the eventual outcome. 

For example, if one is on the "borderline" with regard to encountering 
eventual surge or rotating stall, then a change in the initial amplitude 
can alter the system behavior, from tending to a new equilibrium point to 
undergoing a surge cycle. The range of situations in which the initial 
condition had a significant effect on the ultimate result was not 
investigated in any detail. From the few studies that have been carried 
out, it appears that if B differs by more than about 10 percent from the 
transition value, then the influence of the initial condition is 
substantially diminished. It is felt, however, that the possible 
influence of initial conditions, and thus the question of how to specify 
initial conditions properly, (as well as the details of the nature of 
disturbances at different stages of the surge cycle), is one that must be 
investigated further. 

It is also to be noted that one can always produce changes in system 
behavior by making sufficient changes of initial conditions. For example, 
if one decreases the initial amplitude to 0.00005 (0.01% of mean 
velocity) , the growth of the rotating-stall cell will be so slow during 
the transient that the compressor behavior will be virtually on the 
axisymmetric curve throughout, and a surge will finally result. As 
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discussed previously, however, it is felt that this represents an 
unrealistically weak initial disturbance. We therefore feel somewhat 
encouraged that the calculations will not be very sensitive to the initial 
conditions, although the matter clearly deserves more study. 

Effect of Compressor Length to Radius Ratio 

The other parameter investigated was £ c , the compressor length-to- 
radius ratio. Results are shown in Fig. 7.7, in which trajectories for 
values of % c of 4.0, 6.0, and 8.0 (the base case) are shown. The value 
of B is 1.0 and the initial condition is kept at A(0) = 0.005. As in 
previous figures, the steady-state value of rotating stall amplitude and 
the equilibrium value are also indicated. 

It can be seen that the curves for £ c = 8.0 and 6.0, while 
quantitatively different, have the same qualitative features, namely a 
rapid rise in amplitude and then an essentially quasi-steady tracking 
along the steady-state amplitude curve until the equilibrium value is 
reached. A completely different situation exists for Z c = 4.0, however. 
There is a large discrepancy between the calculated unsteady trajectory 
and the steady-state curve and in fact, the system does not settle down to 
a new equilibrium point. Starting from the initial conditions, the mass 
flow decreases then increases again and the stall cell amplitude decays to 
a negligible value. In other words, the compressor goes into rotating 
stall and then comes out again, returning to the unstalled portion of the 
compressor characteristic. This behavior would thus be a kind of surge 
excursion, rather than an approach to a new system equilibrium in rotating 


stall 
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These calculations show that, in addition to the B parameter, the 
length to diameter ratio, & c , also can affect whether a given system will 
exhibit surge or rotating stall. The reason is that, for the same value of 
B, a decrease in length-to-radius ratio means that the time for the 
formation of a rotating stall cell is relatively longer in proportion to 
the time for a mass flow excursion. This idea appeared in the discussion 
in Section 1 concerning the ratios of different times; it was pointed out 
that the ratio of both the Helmholtz resonator time and the flow-change 
time to the rotating-stall formation time will scale with length-to- 
radius ratio. 

It is relevant here to briefly relate these calculations to previous 
studies. In [5] the only effect that was investigated was the influence 
of the B parameter. In the published results, the "time constant" that 
characterized the rotating-stall formation time was kept constant at two 
rotor revolutions, since the available data was insufficient to determine 
it more precisely. It was, however, stated that some calculations had 
been carried out using different values of that time constant and the 
results for the critical value of B were somewhat different, although over 
the range examined the difference was not great. In the present study, 
there is no need to assume such a time constant, because we are now 
actually calculating the manner in which the compressor output transitions 
from axisymmetric to fully developed rotating stall. 

Comparison with Data on Instantaneous Compressor Performance 

In addition to examining the impacts of B parameter, initial condi- 
tions, and compressor length-to-radius ratio, we have also used the 
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present calculation procedure to help in the interpretation of some 
transient compressor data. The situation is as follows: We have taken 

data from a surge cycle , in the form of overall atmospheric to plenum 
pressure rise versus mass flow, and from this we have attempted to estimate 
the instantaneous compressor pumping characteristic. 

The specific procedure we have followed is to subtract from this 
overall pressure rise a one-dimensional (i.e., axisymmetric) correction 
for the inertial forces due to fluid accelerations during the surge cycle. 
The explicit expression for the correction is: 

*c ' ’ + *c < 7 - 2 > 

Note that this is only valid for axisymmetric flow. The terms given by 
the data are <p and i c d$/d£, and \|> c is calculated. The details of this 
procedure are given in [7,8]. 

In applying this procedure to the data, however, it is also necessary 
to examine the time— resolved data to determine whether the surge includes 
only axisymmetric flow. It was found that the data showed rotating stall, 
not fully developed, toward the end of each period of rapid flow change. 
Although the data are not sufficient to resolve all the details of the 
flow, one must conclude that axisymmetric flow cannot be assumed for the 
entire surge transient that the data represents. 

A sample of the data is shown in Pig. 7.8, with the inertially correc- 
ted characteristics determined by Eq. (7.2). Curves 1 and 2 are found from 
correcting the accelerating and decelerating portions of the surge cycle, 
respectively. Also shown in the figure, by the dot-dash line, is the 
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measured steady-state data for both unstalled and rotating stall operation. 
It can be seen that even though there is some scatter in the calculated 
^ c> due presumably to the differentiation of the experimental data, the 
points fall into two categories: For decelerations, the corrected curve 

lies substantially above the steady-state rotating-stall curve, whereas for 
accelerations (at least for $ > 0.4), the derived pumping characteristic is 
considerably below the steady-state curve. In view of this behavior, the 
question arises whether a single-valued axi symmetric curve is indeed a 
valid representation, or whether a more comlex axisymmetric behavior must 
be invoked. 

To try to understand the situation in more detail, calculations have 
been carried out using the parameters that characterize the transient data. 
These are B = 1.58, l c = 8.0, and a throttle setting that corresponds to a 
throttle curve of K T = 10. The calculations were done both for accelerating 
and decelerating flows. For the decelerations, the initial conditions were: 
$ = 0.5, Y = 0.66, and A(0)W = 0.005. For the accelerating part of the 
cycle the initial conditions were $ = 0.0, ¥ - 0.3, and A(0)W = 0.005. 

The results are shown in Fig. 7.9. The estimated axisymmetric 
characteristic is the solid line and the steady-state rotating stall curve 
is the dashed line. The computed instantaneous compressor performance 
(which corresponds to the quantity shown dashed in Fig. 7.8), is indicated 
by the heavy dash-dot lines. The arrows show the direction of the motion 
of the operating point. 

It can be seen that the computed and the experimental curves have 
strong similarities, at least qualitatively. In particular, the tendency 
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for the compressor pumping characteristic to be above the axisymmetric 
performance in the reverse flow region and below it at high forward flows 
is clearly manifested. 

The reason for this is found from consideration of the growth and 
decay of the rotating stall. During the latter part of the deceleration, 
the instantaneous rotating stall amplitude is much larger than the steady- 
state value (which is in fact zero for negative values of $)• Thus the 
compressor output is closer to that which is characterized by these large 
values of A. Conversely, during the accelerations, for the flows greater 
than $ = 0.4, say, the rotating stall amplitude is again much larger than 
that in steady-state operation, and the performance (i.e., pressure rise) 
also corresponds to that associated with large values of stall cell 
amplitude. 

The main points in the foregoing comparison of experimental data with 
the present computations are thus: 

1 ) The behaviors of the two are qualitatively similar, and 

2) The departure of the compressor pumping characteristic from both 
axisymmetric and steady-state performance is due to the fact that the 
rotating stall amplitude differs considerably from the steady-state 
value. In particular, in the latter part of both deceleration and 
acceleration processes, the stall cell is much "bigger" than one would 
presume from steady-state considerations, and the performance thus 
departs considerably from that of the steady-state. 
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8. DISCUSSION OF THE PRESENT MODEL AND SUGGESTIONS FOR FUTURE WORK 

Although the numerical calculations and the analysis have pointed out 
certain features of the problem which appear to be important, it is clear 
that the present treatment has only made a start on this very complex 
problem. It is thus appropriate to give some discussion about those 
aspects of the actual flow that are included in the theory, those aspects 
that are potentially important but have not been included, as well as some 
of the areas that have been pointed up as needing further investigation. 

We can start by listing those aspects of the unsteady compressor and 
compression system behavior that are accounted for: 

1 ) System dynamic characteristics are included. They are analyzed in a 
lumped parameter fashion. For the situations considered in the present 
report (low Mach number and low frequency) this is a very good approxima- 
tion. As either frequency or Mach number is increased, one may have to 
adopt a more complex description of these components. In particular, there 
is, for high pressure ratio compressors, the possibility of mass storage 
within the compressor itself so that a simple representation of the 
inertance of the flow in the compressor duct may not suffice. In this 
regard, however, it is of interest to note that Mani [10] has achieved some 
quite reasonable comparisons with experiments using very simple represen- 
tations of the system parameters. 

2) Inlet flow field - A description of the upstream flow is included in 
the theory. The flow upstream of the compressor is modeled as two- 
dimensional and unsteady, but irrotational. The upstream annulus is taken 
as having constant area. This assumption can readily be relaxed, if 
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necessary. It is not clear how valid the irrotationality assumption is, 
and this is a subject for investigation. Data is scarce on this point, 
although the comments of Day do give credence to this assumption. 

Perhaps more serious is the assumption that axial and circumferential 
disturbances are related as they would be if they were harmonic. This 
assumption requires study and evaluation. Also, it has been assumed that 
entrance to the inlet guide vane is accomplished without local head loss. 

3) Exit flow field - A description of the exit flow field is included in 
the theory. A more important assumption is that of linearity in the exit 
pressure perturbations. It is clear that this cannot be correct, in 
general, although the degree to which it is in error and, more importantly, 
its effect on the predictions of the theory is not clear. Previous 
theories of steady-state rotating stall, in particular [3] and [11] which 
have made this assumption, have had good success with predicting at least 
some of the features of the stall cell flow field. 

4) The model of the throttle characteristic that is included is quite 
adequate for low pressure ratios. For higher pressure ratio systems, 
there will be no problem with introducing a compressible version of this, 
as was done by Wenzel [12]. 

5) Throttle transients (e.g., ramp closures) have not been included. 
However there should be no difficulty in applying the theory to account 
for these. This is of interest since it has been found that throttle 
transients can have a significant effect on system behavior. 

6) System hysteresis is included in the theory. This appears naturally 
since it is set by the stability of the system at the various intersection 
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points of throttle and compressor curves. When closing the throttle (into 
stall) the compressor curve is the axisymmetric one. When opening the 
throttle from operation in fully developed rotating stall, the compressor 
curve is the rotating stall one, and there can thus be a difference in the 
throttle settings at which transition from unstall to stall and from stall 
to unstall occur. 

7) Compressor geometry is represented by the axisymmetric characteristic. 

8) Unsteady blade row response (internal aerodynamic lag of the compressor) 
is included in the theory. It is done in a rudimentary manner, because at 
present one does not know how much detail is needed to provide a satisfactory 
description. 

From the above it does appear that all of the elements that one might 
think necessary for an analytical description of general post-stall tran- 
sients have been included. Having said this, however, it should be 
emphasized that their inclusion has, in some instances, been done in a 
quite approximate manner, with an impact on the overall prediction which is 
not well understood. In this connection, we can refer back to the basic 
analysis [3] from which this study begins. The equations developed there 
will always result in axial velocity profiles that are like those shown in 
Fig. 3.4, for example. They will not yield the 11 top hat" type of profile 
that seems to be more characteristic of compressors that operate in rotating 
stall. At present, the reason for this discrepancy is unclear, but the 
fact that it exists suggests that there are aspects of the theory which 
should be improved, even for low speed, low pressure-ratio situations. 

We have described those features that are included, and we now turn 



68 


to those that have not been, but which we think will be important to 
include for problems of practical interest. 

1 ) Effects of compressibility on compressor performance and inlet and 
exit flow fields. 

2) Front-to-rear mismatching effects due to off-design operation in the 
compressor. 

3) Heat addition in the burner. 

4) Effects of inlet flow distortion on post-stall transients. 

5) Engine matching effects (non-constant speed, for example). 

6) Other compression system components. 

7) Steady-state hysteresis in the axisymmetric compressor characteristic, 
if such exists. 

Research is necessary to elucidate the influence of the effects named 
above. In addition to these, however, there are other topics that must be 
explored if one is to improve the quantitative capabilities of the theory. 
The first of these is the development of more exact methods for the 
calculations. The numerical results given in the present report are based 
on the one-term Galerkin procedure; as mentioned, this was all that could 
be carried out in the limited time available. In particular, the type of 
phenomena that are being investigated would appear to be very well-suited 
for use of spectral method calculations. Such procedures would be able to 
directly use any Fourier series representations of the inlet and exit flow 
fields that are developed. At present, we are examining ways to carry out 
this calculation, and so we will only mention that this is one item that 


should be addressed. 
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Another item is concerned with the rotating stall performance. The 
question of the behavior of the compressor at or near the point of 
instability is by no means settled. The present theoretical analysis 
predicts that the flow in the compressor will become unstable at the peak 
of the steady-state pressure rise curve. In practice, compressors are 
observed to become unstable on the negatively sloped part of the curve. 

While interstage (one-dimensional), volume effects have been invoked to 
account for this in high pressure ratio compressors, this is observed at 
low pressure ratios also, and here the interstage effects are very small. 
Research is needed on this topic, both of an experimental as well as a 
theoretical nature. 

The effect of inlet conditions has been mentioned before, and the 
point will not be belabored here. We will merely note that this is a new 
aspect of the problem, and one on which there is little data. 

Effects of initial disturbance features on the subsequent transients 
have received only cursory attention in this report. Only a rotating-s tall- 
like harmonic disturbance was considered, and it began when the throttle 
was set and held precisely at the neutral point. Obviously, a myriad of 
other possibilities need study; disturbances which are statistically 
defined fluctuations, or which arise from inlet distortions of various 
types, or which arise from asymmetric burner pulses, are all possibilities 
for transient initiation, perhaps in combination with symmetric, or surge- 
like, pulses arising either naturally or from control actions. Research on 
these questions of disturbance initiation should be helpful to guide 


experimental tests 
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9. SUMMARY AND CONCLUSIONS 

1 ) An approximate theory has been developed for post-stall transients in 
axial compression systems. 

2) The analysis includes a two-dimensional unsteady representation of the 
compressor flow field, together with a description of the overall dynamic 
response of the system. A system of three simultaneous first-order 
ordinary differential equations are derived, through which surge-like and 
rotating-s tall-like oscillations are coupled. 

3) Examination of the coupled equations shows that surge and rotating 
stall can each exist in "pure " or equilibrium form; however, rotating 
stall cannot evolve without inducing some surge-like unsteadiness in the 
process. During either equilibrium surge or rotating stall, disturbances 
of the other family will tend to die out; that is, "pure" modes tend to be 
stable • 

4) Transients are examined, resulting from a specified initial 
disturbance. They all have unsteady mass flow through the compressor, 
thus including situations in which rotating stall is growing or decaying. 
This growth or decay is calculated, along with the instantaneous system 
operating point. To the authors* knowledge, this is the first time that 
this has been done. 

5) The instantaneous rotating stall cell amplitude is found to have a 
significant effect on the instantaneous compressor pumping characteristic. 
This, in turn, affects the overall system response. 

6) The theory extends previous analysis about the effect of the B para- 
meter (B « (U/2a s )/V/AL c ) on post-stall transients, and adds details 
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of the transient process; how rotating stall decays for large B, and how it 
settles to a fixed level when B is less than the critical value. 

7) Other parameters besides B, in particular the compressor length to 
radius ratio x, c = L/R, can also have a strong effect on the system 
response. A limited study shows that (for a given value of B) compressors 
of shorter length to radius ratio are more likely to exhibit surge than 
rotating stall. The initial conditions concerning stall cell amplitude at 
the start of the transient are presumably also of importance, and should 
carefully be studied. 

8) The rotating-stall cell amplitude during unsteady flow is different 
from that during steady-state operation in rotating stall. Consequently, 
the instantaneous compressor performance during a system transient can 
differ considerably from the characteristic measured during steady-state 
rotating stall. 

9) The numerical results that are presented show qualitative agreement 
with the existing (but scarce) experimental data concerning the nature of 
the flow field during this type of transient. 

10) Based on the results of the analysis, recommendations are made 
concerning future work in the area of stagnation stall. 
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FIGURES 



Figure 2.1 Schematic of compressor/compression system geometry. 
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Figure 3.1 Notation used for compressor performance characteristics. 

Total-to-static pressure rise (\Jj c ) vs. flow coefficient (<f>) 
in absence of rotating stall. Operating point in rotating 
stall is <5 above ty c . 
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Figure 5.1 Schematic of compression system showing nondimensional lengths. 
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Figure 6.2 Amplitude of rotating stall during a surge cycle, as indicated 
by terms in Eq. (6.9), with 6=0. Different choices of £ A 
provide different starting points. Net amplification occurs 
in shaded zone. 



Figure 6.3 Sketch of surge cycle showing regions in which growth of 

rotating stall can occur. Numbered points correspond to lines 
of Figure 6.2. 
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Figure 7.2 Transient compression system response (t vs. $) for different 
values of B, (£ = 8.0/ A(0)W = 0.005). 
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Figure 7.4 Evolution of rotating stall amplitude during system transient - effect of initial 
conditions; stall cell ~ A(£)Wsin(0-f£) , (£ = 8.0, B = 1.0). 
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Figure 7,6 Instantaneous compressor pumping characteristic during system 
transient - effect of initial conditions (£ c = 8.0, B = 1.0). 
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Figure 7.7 Evolution of rotating stall amplitude during system transient - effect of length- to-radius ratio; 
stall cell~A(£)Wsin(0-f£) , (A(O)W = 0.005, B = 1.0). 
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Figure 7.8 Instantaneous compressor pumping characteristic derived from 
surge cycle data; measured steady-state rotating stall curve 
also shown for reference. 
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Figure 7.9 Calculated instantaneous compressor pumping characteristic 
(B = 1.58, i = 8.0, A (0) W = 0.005). 






Figure B.1 Calculations of 6 for numerical (dot— dash) , one— harmonic 

Galerkin (dash), and two— harmonic Galerkin (short dash for 
shallow, solid for steep) diagrams. I, II, III denote three 
solutions for latter case. 
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APPENDIX A 

The "dh/dQ = -g 11 Approximation 

The relation between h and g is discussed in Ref. [3] and more fully 
in Appendix A of Ref. [4]. No simple exact relationships apply except in 
special circumstances. When the entrance duct is infinite in length and of 
constant width, an integral transformation connects g and h: 


h = 


TT 


2tt 

g (z)ln 

0 





(A . 1 ) 


Alternatively, since g(6), being periodic, has a Fourier series of the form 


g = I (a sin n0 + b cos n0 ) (A. 2) 

n=1 n n 


The corresponding series for h would be 


h = I (a cos n0 
n=1 n 


b sin n9 ) 
n 


(A. 3) 


Eqs. (A. 1 ) and (A.2-A.3) are equivalent. The average of g ( e ) over a cycle 
vanishes by definition. Therefore, there is no n = 0 term in Eq. (A.2). 

In Ref. [3] part II, it is shown that an exact solution (i.e., a 
particular set of a n and b n ) applies when G(g) is purely parabolic, 
with zero slope at the throttle point g = 0. The solution may be written 

i0 


h + ig = K 


1 + rje 


ie 


(A. 4) 


A 

where n is an amplitude parameter. 

If rotating stall oscillations have a small amplitude (Ref. [3] part 
I), then only the first terms of Eqs. (A.2, A. 3) are needed. An obvious 
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example is provided by Eq. (A. 4), where n may be set equal to 0 for 
small amplitude. In that case, 

dh/de « -g(0) (A. 5) 

would be true. In Ref. [3] part III , Eq. (A. 5) was used in Eq. (3.7) for 
limit-cycle calculations, even though it is strictly true only for weak 
disturbances from uniform flow. 

In Appendix A of Ref. [4] , it was shown that an analysis based on Eq. 
(A. 5) may be validated for two harmonics of Eqs. (A. 2, A. 3): One defines 

a new function z(0) as follows: 

g = -dz/de ; h = Az - Bd 2 z/d0 2 (A. 6) 

Thus, from Eq. (A. 2), 

oo 1 

z(0) = £ — (a cos n0 - b sin n0 ) 

n=1 n n n 


and hence, Eq. (A. 6) implies that 


h = E ( — + nB ) (a cos no - b sin n0 ) 
n«1 n n n 


Comparing Eqs. (A. 7) and (A. 3), one sees that A and B, until now 
arbitrary, may be chosen to make two harmonics agree (that is, make the 
first parenthesis equal to 1 for two different values of n). If, for 
example, the first two harmonics are to agree, then it is easily seen that 



are the proper choices. 

When Eq. (A. 6) is introduced into Eq. (3.7), one finds 
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~ 2 2 ^ 

(X - Bmf )d z/d0 + mfAz + G(-dz/d0) -6=0 (A. 9) 

This is exactly the equation that would apply (for h instead of z) if Eq. 
(A. 5) were used, provided X and f are redefined. Thus, one may conclude 
that the approximate relation dh/d0 = -g, though not exact, can be adopted 
for nonlinear calculations with some confidence, and we will do so in this 
report. We will have in mind that the procedure described by Eq. (A. 6) 
would improve accuracy, but we leave that refinement for future studies. 
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APPENDIX B 

Two-Harmonic Galerkin Method 

In general, the Galerkin method is based on some suitable sequence of 
functions which, after the corresponding parameters are determined, can 
well represent the wave form of the oscillation. If the sequence is taken 
to be harmonic, of the following form, in which amplitudes have been 
scaled with the diagram semi-width W: 

g = WA* sine + WB* sin2(e-u) (B.l) 

then the Galerkin procedure is identical to a method of straightforward 
solution by Fourier series. Here, we will carry out the harmonic analysis 
for pure rotating stall to second order, showing how the three unknowns A*, 
B*, and u are found; the extension to higher orders will be obvious. 

Using the dh/de = -g approximation, we find 

h = WA* cose + 1/2 WB* cos2(0-u) (B.2) 

The requirements of Eqs. (3.10) and (3.11) are automatically met by use of 
Eqs. (B.l) and (B.2). We note that if Eq. (A. 3) were applicable instead 
of Eq. (A. 5), then Eq. (B.2) would be replaced by 

h = WA* cose + WB* cos2(e-u) (B.3) 

It remains to satisfy the differential equation, Eq. (3.9), with G 
specified by Eq. (3.16). Following the Galerkin procedure, we substitute 
Eq. (B.l), and, denoting the left side of Eq. (3.9) as L(e), we find 
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^L(9) = (-X + mf)A*sin9 + (— 4X + mf )B*sin2( 9-u) + ... 

2 

+ H (1 ■ ^ " 2 W [A * sine + B*sin2(9-u) + ...] 

W 

- [A*sln9 + B*sin2(9-u) + . ,.] 2 

x [A*cos9 + B*cos2 ( 9-u ) + ...] (B.4) 


Then, we multiply L ( 9 ) by the representational functions sin9 , cos9 , 
sin29 , cos29, . . ., and require that the integrals of each of those 
products vanish over a cycle. To second order, we have four equations for 
the unknowns X, A*, B*, and u. With each step in order, two new moment 
equations and two new unknowns (amplitude and phase angle) are added. 

Using the cubic form Eq. (3.16), and recalling that the average of g 
vanishes, the performance effect from Eq. (3.12) is 


6 


_J 

4-itW^ 


[36 


2ir 2 

g d9 + 

0 


g 3 dg ] 

0 


(B.5) 


into which the profile Eq. (B.1) may be substituted. By Eq. (3.5), the 
total in-stall performance y is i)> c ($) + 6. Eq. (3.13) with $ = W + 6 
gives 


Y = 6 + i|i + H(1 
c 

o 


3 B 11 

2 W “ 2 3 

W 


(B.6) 


The two-harmonic analysis has been carried out, with B*, X, A*, and \y 
found numerically from a set of nonlinear algebraic equations. This was 
done for a variety of values of a, H/W, and B- In addition, the effect of 
the dh/de* = -g approximation was investigated by keeping or not keeping 
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the factor 1/2 in Eq. (B.2). 

Two helpful conclusions can be drawn from the results, namely that the 
wave profile g(0) and performance effect 6 are only weakly affected by the 
lag parameter a, and that the dh/d0 1 = -g approximation is quite accurate. 
These results tend to support the one-harmonic analysis, in which those two 
features play no role at all. 

The diagram steepness (H/W) is indicated to be a very important 
parameter, however, even though it too is absent from the one-harmonic 
result {except as a trivial matter of scale). The effect of H/W is best 
illustrated in terms of performance effect (6), shown on Fig. B.l. There, 
the right side of a cubic diagram ( ) is illustrated for H/W = 0.67 and 
0.13, and below them are shown certain results for 6. As we already know, 
the one-harmonic result is the parabola, shown dashed, whatever the value 
of H/W. The numerical result (for H/W = 0.72) shown dot-dashed has a 
different shape, with a much more abrupt change near recovery. This is 
presumably dependent on H/W. 

The two-harmonic solution for H/W = 0.13 is shown as the short-dashed 
line. Even for a "shallow" diagram, there is a shape change of the 
function 6(0). Of course, as H/W approaches zero, the one-harmonic result 
is approached as a limit. As H/W is increased beyond about 0.2, the 
two-harmonic result is unfortunately not unique. Of course, one does not 
know that the exact solution is in fact unique; however, the disjointed 
character of the three solutions (indicated as solid lines) for H/W = 0.67 
suggests that the harmonic Galerkin method cannot be stopped at second 
order. Many harmonics would presumably be needed to achieve agreement 
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with the exact solution, whether or not the solution is truly unique. 

Further attempts to find an improvement over the one-harmonic Galerkin 
analysis which is still relatively simple have failed. If it were purely 
a matter of representing the numerical profiles (Figs. 3.4a-f), then the 
method of matched asymptotic expansion by which Cole (Ref. [13]) treated 
the van der Pol problem could be effectively used. However, that method 
in effect assumes rapid changes ("relaxation") in parts of the wave 
process, and we cannot accept that feature if we wish to consider the 
initial growth of a rotating-stall disturbance; in that case, a single 
harmonic wave would correctly describe the first appearance of an 
oscillation. It seems that improvements beyond the one-harmonic Galerkin 
method will have to be numerical. 
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NOTATION 

A amplitude function of first-harmonic angular disturbance, Eq. (5.50) 

A* amplitude of first-harmonic angular disturbance in pure rotating stall, 
Eq. (3.17) 

A 

A disturbance parameter, Eq, (5.33) 

A c compressor duct area 

a reciprocal time-lag parameter of blade passage, Eq. (3.4) 
a n Fourier coefficient of axial flow disturbance 

a s sound speed 

B (U/2a s ) yv/AL c 

B* amplitude of second-harmonic angular disturbance in pure rotating 
stall, Eq. (B .2 ) 

A 

B disturbance parameter, Eq. (5.33) 

b blade chord 

b n Fourier coefficient of axial flow disturbance 

C x axial flow velocity 

D mean wheel diameter 

E number of order 1, Eq. (1.5) 

F pressure-rise coefficient in blade passage, Eq. (2.4) 

F t throttle characteristic function, Eq. (5.36) 

F^ 1 inverse of throttle function 

f speed coefficient, relative to laboratory, of propagation of 
angular disturbance 

f Q value of f for harmonic-wave rotating stall, Eq. (5.55) 

G compressor characteristic disturbance function, Eq. (3.6) 
g disturbance of axial flow coefficient 

H semi-height of cubic axisymmetric characteristic, Fig. 3.2 
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h circumferential velocity coefficient 

J square of amplitude of angular disturbance of axial-flow coefficient, 
Eq. (5.57) 

J Q initial value of J 

J e value of J for fully-developed rotating stall at the existing 

average axial-flow coefficient, Eq. (6.1) 

K coefficient, Eq. (1.2) 

A 

K parameter of exact solution, Eq. (A. 4) 

Kq loss coefficient at IGV entrance, Eq. (5.6) 

K T parabolic throttle coefficient, Eq. (4.3) 
k factor accounting for interrow spacing, Eq. (5.21) 

k T linear throttle coefficient, Eq. (4.4) 

L left side of an equation; zero exactly, made to vanish on certain 
averages in a Galerkin method 

L c total effective length of compressor and ducts, Eq. (1.2) 

L r axial length of a blade row, Eq. (5.21) 

«. c total length of compressor and ducts, in wheel radii, Eq. (1.7) 

Jtj length of entrance duct, in wheel radii, Fig. 5.1 
length of exit duct, in wheel radii, Fig. 5.1 
«, T length of throttle duct, in wheel radii, Fig. 5.1 
m duct-flow parameter of rotating stall, Eq. (3.2) 

N number of stages of core compressor 
P pressure coefficient, Eq. (5.13) 

P 0 static pressure at entrance to IGV, Fig. 2.1 

Pi static pressure at entrance to core compressor. Fig. 2.1 

P E static pressure at compressor exit, Fig. 2.1 

P S static pressure at end of exit (diffuser) duct, and pressure in 

the plenum. Fig. 2.1 
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P T total pressure ahead of entrance and following the throttle duct, 

Fig. 2.1 

r time-dependent phase angle, Eq. (5.50) 
t time 

U wheel speed at mean diameter 

u velocity in circumferential direction 

Vp volume of plenum, Fig. 5.1 

v velocity in axial direction 

W semi-width of cubic characteristic, Fig. 3.2 

Y disturbance potential at compressor entrance, Eq. (5.27) 

A 

Y function defined in Eq. (5.33) 

y axial distance, Fig. 2.1 

Z number of order 1, Eq. (1.3) 
z function defined in Eq. (B.3) 

B location of throttle setting, ($/W)-1 

Y stagger angle of blades, 50% reaction 
AP compressor pressure rise, Eq. (1.2) 

6 performance effect of pure rotating stall, Eq. (3.5) 
n axial distance measured in wheel radii, Eq. (2.2) 

t\ parameter of exact solution, Eq. (A. 4) 

0 angular coordinate around wheel, Fig. 2.1 

0* angular coordinate around wheel, following disturbance, Eq. (5.45) 

X rotating-stall parameter, Eq. (3.2) 

u phase angle defined in Eq. (B.1) 

£ time, referred to time for wheel to rotate one radian 

£ e dimensionless time for development of equilibrium rotating stall, Eq. (6.4) 
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5a defines time of start of angular disturbance, Eq. (6.8) 
p density 

t coefficient of pressure-rise lag, Eq. (2.4) 

t c flow change time scale, Eq. (1.2) 

r e plenum emptying time, Eq. (1.5) 

th Helmholtz resonator time, Eq. (1.4) 

tr rotating-stall formation time scale, Eq. (1.1) 

$ axial flow coefficient, averaged over angle; velocity divided 
by wheel speed 

$ flow coefficient averaged over both angle and time 

$fji flow coefficient of throttle duct, referred to entrance-duct 
area, Eq. (5.35) 

<|> axial flow coefficient C x /U; variable in angle and time 
<j> velocity potential in entrance duct, Eq. (5.7) 

$' disturbance velocity potential, Eq. (5.10) 

Y total-to-static pressure-rise coefficient, Eq. (3.2) 

i|> pressure-rise coefficient, AP/pU 2 

axisymmetric pressure-rise coefficient, Eq. (3.2) 

♦ co shut-off value of axisymmetric characteristic, Fig. 3.2 
Ur Helmholtz frequency, Eq. (1.4) 

Subscripts 

0 at the entrance to the compressor 

E at the exit of the compressor 
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